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I.   INTRODUCTION 

A  well-known  programming  problem  is  the  traveling- salesman 
problem.   As  classically  stated,  a  salesman  must  start  from  his  home- 
town, travel  to  each  of  n-1  other  cities,  and  return  home.   He  is  re- 
quired to  enter  each  of  the  n  cities  exactly  once  and  leave  each  city 
exactly  once,  with  the  understanding  that  his  entire  trip  is  one 
continuous  circuit.   There  is  no  requirement  that  the  distances  be- 
tween cities  be  symmetric,  that  is,  the  distance  from  city  i  to  city  j 
need  not  be  the  same  as  that  from  j  to  i.   The  salesman's  objective  is 
to  select  a  sequence  of  cities  in  such  a  way  as  to  minimize  the  total 
distance  traveled. 

There  are  several  algorithms  which  solve  the  traveling- salesman 
problem,  but  none  offers  the  last  word  in  terms  of  speed  and  ease  of 

computation.   The  purpose  of  this  investigation  was  to  examine  a  newly 

1 
proposed  solution  procedure  developed  by  Professor  Harold  Greenberg. 

His  algorithm  looked  very  promising  in  terms  of  hand  computation  on 

small  problems,  but  this  study  focused  on  its  desirability  for 

handling  larger  problems,  such  as  the  forty-city  case.   The  greatest 

effort  went  into  programming  the  algorithm  for  the  computer,  then 

letting  the  computer  serve  as  the  primary  tool  of  analysis. 


1Greenberg,  Harold,  Professor,  Department  of  Operations  Analysis, 
Naval  Postgraduate  School,  Monterey,  California 


A.   DEFINITIONS 

The  distance  (cost)  between  city  i  and  city  j  is  denoted  c^  for 

i  =  1, . . . ,n  and  j  =  l,...n.   In  general,  c-.  need  not  equal  c...   The 

-J  J  -*- 

decision  variable  x. .  equals  one  if  the  salesman  is  to  travel  directly 
from  city  i  to  city  j  and  equals  zero,  otherwise.   Since  traveling 
directly  from  city  i  to  city  i  is  not  meaningful,  x. .  must  always 
equal  zero  for  all  i.   Most  solution  procedures  meet  this  requirement 
by  assigning  arbitrarily  large  values  to  c.  for  all  i, 

Associated  with  each  city  is  the  number  of  steps  required  by  a 
solution  to  arrive  at  that  city.   For  example,  if  the  salesman  is  to 
go  from  city  A  to  B  to  C  and  back  to  A,  he  arrives  at  B  in  one  step,  C 
in  two  steps,  and  back  to  A  in  three  steps.   Unless  otherwise  stated, 
the  following  discussion  takes  city  one  to  be  the  city  of  origin  and 
city  n  to  be  the  last  city  visited  prior  to  returning  to  city  one. 

A  solution  must  provide  a  tour  in  order  to  be  a  traveling- salesman 
solution.   A  tour  is  defined  as  a  sequence  of  transitions  from  city 
one  through  all  other  cities  and  back  to  city  one,  such  that  each  city 
is  entered  exactly  once  and  left  exactly  once,  and  it  must  take  exactly 
n  steps  to  return  to  city  one.   The  latter  clause  guarantees  there  are 
no  discontinuities  in  a  tour.   Deleting  this  clause  yields  the  defini- 
tion of  a  subtour.   It  is  possible  for  a  candidate  solution  to  contain 
as  many  as  n/2  subtours  and  still  meet  the  requirement  that  each  city 
is  entered  and  left  exactly  once.   To  illustrate  the  concept  of  tours 
and  subtours,  consider  these  solutions  to  the  four  city  problem:   city 
1  to  3  to  4  to  2  to  1  provides  a  tour;  whereas  city  1  to  2  to  1  and 
city  3  to  4  to  3  provides  two  subtours. 


B.   ALTERNATE  FORMULATIONS 

The  traveling-salesman  problem  is  essentially  a  combinatorial 
problem  and  there  is  no  one  way  to  solve  it.   The  most  direct  approach 
is  to  select  a  starting  point,  evaluate  the  distance  (cost)  over  each 
of  the  (n-1)!  possible  tours,  and  then  choose  the  tour(s)  with  minimum 
total  cost.   A  dynamic-programming  formulation  reduces  the  number  of 
tours  one  must  consider  but  the  tables  required  can  quickly  exceed 
computer  storage  capacity.   Branch-and-bound  techniques  can  also  be 
employed  but  suffer  from  similar  limitations  as  dynamic  programming. 
A  fourth  possibility  might  be  to  construct  a  solution  procedure  using 
an  algorithm  for  finding  minimal-cost  flows  in  a  network. 

The  traveling- salesman  problem  can  be  stated  as  an  integer  program- 
ing problem.   Hadley  /~Ref .  6.7  and  Dantzig  /~Ref.  17  present  essential- 
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ly  the  same  formulation.   Dantzig  lets  x.,   =  1  or  0  according  to 

th 
whether  or  not  the  salesman  travels  from  city  i  to  city  j  on  the  t 

step,  where  i,j,t  =  l,...,n.   Then,  defining  x.  .   ,,  =  x..,   he 

i,J5n+l    ijl, 

states  the  problems  as  follows: 

minimize  ^=*"'   ex    =  z 

i,j,t   ij   ijt 

subject  to  ^-r  x    =  ,   x         (jjt=  l,...,n) 

i   ijt    k   j,k,t+l 

and  ^>_  x      =  1   (i  =  l,...,n) 

j»t   i,j,t 

and  x. .   =  0  or  1  for  all  i,i,t. 

ijt 

The  first  set  of  constraints  requires  that  if  the  salesman  arrives  at 
city  j  in  t   steps,  he  must  leave  that  city  on  step  t+1.   The  second 

sets  requires  that  the  salesman  leave  each  city  i  once  and  only  once. 

3  2 

Solving  this  system  of  n  +1  variables  and  n  +  n  equations  with  an 

integer  programming  algorithm  yields  a  minimum-distance  tour. 


The  integer  programming  approach  has  several  shortcomings.   The 
third  subscript  t  on  x  greatly  increases  the  number  of  variables,  and 
the  number  of  constraints  in  the  formulation  is  very  large.   Further- 
more, integer  programming  algorithms  normally  require  many  more  itera- 
tions than  ordinary  linear  programming  techniques „ 

There  is  a  general  class  of  indirect  solution  procedures  that 

takes  advantage  of  a  special  feature  of  the  traveling-salesman  problem. 

If  one  temporarily  ignores  the  fact  that  a  solution  must  provide  a 

tour,  the  salesman's  task  can  be  formulated  as  a  simple  assignment 

problem,  perhaps  the  most  readily  soluble  of  all  linear  programs. 

The  assignment  formulation  is  achieved  by  letting  x   =  1  or  0 

ij 

according  to  whether  or  not  the  salesman  is  to  travel  from  city  i  to 

i.   Let  c   be  arbitarily  large  for  all  i.   Then, 
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minimize  ^ —  *         c   x 

i    J    ij   ij 


subject  to  ^   x. .  =  1    (j  +  l,,..,n) 

and  *i    x   =1     (i  =  l,...,n) 

J    ij 

2 
Notice  that  there  are  only  n  variables  and  2n  constraints.   One 

constraint  may  be  dropped  because  exactly  2n-l  are  linearly  independent, 
Furthermore,  since  the  assignment  problem  is  a  form  of  the  transporta- 
tion problem,  and  since  all  restraining  expressions  equal  one,  any 
solution  should  yield  0  or  1  values  for  the  x. .'s  /_Ref.  2/ .   Thus, 
there  is  no  need  to  add  the  integer  constraint  to  the  formulation 
above.   Possible  solution  procedures  include  the  simplex  method,  the 


familiar   transportation  tableau,  or  an  algorithm  for  finding  minimal- 
cost  flow  in  a  network  /_Ref.  4/. 

An  assignment  outcome  does  not  necessarily  provide  a  traveling- 
salesman  solution,  but  it  does  satisfy  the  requirements  that  each 
city  is  entered  exactly  once  and  departed  exactly  once  and  that  x. .=  0 
for  all  i.   Clearly,  since  the  salesman's  problem  includes  the  added 
constraint  that  any  set  of  values  for  the  x-.'s  must  provide  a  tour, 
his  minimum-distance  solution  must  be  at  least  as  great  as  the  minimum 
assignment  distance. 

The  so-called  indirect  procedures,  alluded  to  above,  typically  use 
an  assignment  solution  as  a  starting  point.   The  optimal  value  of  the 
assignment  objective  function  provides  a  lower  bound  to  any  traveling- 
salesman  solution.   A  non-optimal  assignment  solution  may  still  provide 
a  lower  bound  while  being  closer  in  magnitude  to  the  salesman's  minimum 
total  distance.   In  the  sense  that  it  is  closer,  the  latter  lower 
bound  is  said  to  be  "better"  than  the  former.   In  practice,  however, 
it  is  difficult  to  find  a  value  larger  than  the  optimal  assignment 
solution  and  prove  that  it  is,  in  fact,  a  lower  bound  to  the  traveling 
salesman  problem. 

Professor  Greenberg's  algorithm  uses  the  optimal  assignment  solu- 
tion as  its  point  of  departure  because  it  provides  a  guaranteed  lower 
bound  and  because  special  use  is  made  of  the  coefficient  vectors  in 
the  final  simplex  tableau.   After  the  solution  of  the  initial  assign- 
ment problem,  the  algorithm  simply  manipulates  the  coefficient 
vectors  until  a  minimum-cost  tour  is  achieved.   This  procedure  is 
explained  in  detail  in  the  next  chapter. 


W.  L.  Eastman  also  begins  with  the  optimal  assignment  solution  and 

if  it  should  provide  a  tour,  then  the  problem  is  solved.   If,  in  the 

more  likely  event,  it  does  not,  he  selects  the  subtour  with  the  least 

number  of  variables  and  sets  each  equal  to  zero,  one  at  a  time.   For 

each  case,  he  resolves  the  original  assignment  problem  but  with  the 

added  constraint.   In  practice  this  can  be  achieved  by  setting 

c. .  =  &&    corresponding  to  the  x. .  required  to  be  zero.   If  a  tour  is 
ij  1J 

not  obtained,  he  branches  from  that  assignment  problem  with  the 
smallest-valued  objective  function.   Again  the  least-membered  subtour 
is  chosen  to  create  new  assignment  problems,  and  so  forth,  until  a 
tour  is  achieved,   /Ref.  8/ 

J.  D.  C.  Little  and  others  have  devised  a  branch-and-bound 
algorithm  that  starts  with  a  lower  bound,  not  to  the  traveling  salesman 
problem,  but  to  the  assignment  problem.   Instead  of  obtaining  the 
optimal  solution  to  the  assignment  problem  at  each  step,  the  smallest 
element  is  subtracted  from  each  row  of  the  cost  matrix,  and  then  the 
smallest  element  from  each  column  of  the  result.   The  sum  of  the  re- 
maining elements  is  said  to  be  a  valid  lower  bound  on  the  optimal 
assignment  solution.   At  each  branching  step  two  new  problems  are 
created  corresponding  to  x,-  •  =0  and  x.  .  =  1.   The  x.  .  chosen  for  each 
branching  step  is  the  one  that  yields  as  large  a  bound  as  possible 
when  set  equal  to  zero.   Branching  proceeds  until  a  tour  is  obtained. 
Little's  method  normally  requires  many  more  branches  than  does  Eastman's, 
but  there  is  less  computation  at  each  step.   Some  performance  results 
are  available  on  Little's  algorithm,,   The  mean  execution  time,  on  an 
IBM  7090  computer  for  100  thirty-city  problems  was  58.5  seconds  and  the 
average  for  five  forty-city  problems  was  8.37  minutes.   /Ref.  9/ 
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C.   PLAN  OF  THE  REPORT 

Chapter  II  presents  a  detailed  description  of  the  solution  procedure 
suggested  by  Professor  Greenberg.   His  proposed  algorithm  is  complete 
in  the  sense  that  the  criteria  for  proceeding  from  one  step  to  the 
next  are  completely  specified  and  guarantee  an  optimal  solution.   In 
some  instances,  however,  it  is  not  readily  apparent  how  one  can  most 
efficiently  meet  these  criteria.   Chapter  III  deals  with  these  ques- 
tions of  technique  as  it  describes  how  the  algorithm  was  programmed 
for  the  computer.   When  there  were  alternate  computational  approaches 
possible,  an  attempt  was  made  to  find  a  good  one.   There  is  no  claim 
that  the  particular  techniques  devised  to  meet  the  algorithm  criteria 
are  the  ones  Professor  Greenberg,  himself,  would  have  chosen. 
Chapter  IV  outlines  the  salient  results  of  computational  experience 
with  the  computer  program.   The  algorithm  and  the  program  are  evaluated 
in  the  light  of  this  experience. 
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II.   THE  ALGORITHM 

An  algorithm  is  defined  by  Webster  as  "a  rule  or  procedure  for 
solving  a  mathematical  problem  that  frequently  involves  repetition  of 
an  operation."   This  chapter  seeks  to  present  Professor  Greenberg's 
algorithm  for  solving  the  traveling-salesman  problem  in  the  logical 
sequence  of  its  development,  concluding  with  a  concise  statement  of 
the  algorithm,  itself. 

A.   THE  TRAVELING- SALESMAN  PROBLEM 

A  traveling  salesman  is  faced  with  the  problem  of  finding  the 
shortest  route  (or  the  least-cost  route)  through  n  cities.   Further- 
more, he  must  meet  the  following  conditions:   1)  he  must  enter  each 
city  once  and  only  once,  2)  he  must  depart  from  each  city  once  and 
only  once,  and  3)  his  route  must  be  a  tour,  that  is,  there  must  be  no 
discontinuities  in  his  route.   There  are  (n-1)!  possible  solutions  to 
this  problem,  since  the  salesman  can  go  to  any  of  n-1  cities  on  his 

first  step,  then  to  any  of  the  remaining  n-2  cities  on  his  second 

th 
step,  and  so  forth,  until  he  returns  to  his  starting  point  on  the  n 

step.   Optimal  solutions,  of  course,  are  found  within  this  set  of 
(n-1) !  tours. 

Clearly,  the  above  is  merely  an  allegorical  statement  of  a  problem 
that  has  relevance  to  a  variety  of  situations.   For  example,  suppose 
fixed  electronic  components  must  be  connected  in  a  continuous  circuit 
such  that  the  least  amount  of  wire  is  required.   Disregard  of  condi- 
tions (1)  and  (2)  might  leave  some  components  out  of  the  circuit 
entirely  or,  on  the  other  hand,  cause  short  circuits.   Violation  of 
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condition  (3)  would  result  in  two  or  more  subcircuits,  with  current 
flowing  in  only  one. 

B.   THE  ASSIGNMENT  FORMULATION 

Modifying  the  traveling-salesman  problem  by  disregarding  the  third 
requirement  above,  yields  a  problem  which  is  readily  adaptable  to  the 
assignment  formulation,  perhaps  the  most  easily  soluble  of  all  program- 
ming problems.   The  assignment  problem  is  usually  stated  in  terms  of 
assigning  n  men  to  n  jobs  in  such  a  way  that  overall  productivity  is 
maximized,  assuming  that  a  productivity  index  can  be  associated  with 
each  man-job  combination.   There  are  n!  possible  solutions  to  this 
problem,  since  the  first  man  can  be  assigned  to  any  of  the  n  jobs,  the 
second  man  to  any  of  the  remaining  n-1  jobs,  and  so  forth. 

Dropping  the  requirement,  then,  that  a  solution  must  provide  a 
tour,  the  modified  salesman  problem  is  to  assign  n  cities  to  n  cities 
(the  same  n  cities)  such  that  the  sum  of  city-to-city  distances  is 
minimized.   The  possibility  of  assigning  city  i  to  itself  is  precluded 
by  giving  each  c. .  an  arbitrarily  large  value.   The  possibility  of 
other  types  of  subtours  appearing  in  a  solution,  however,  cannot  be 
restricted  without  destroying  the  assignment  formulation. 

Since  no  city  is  ever  assigned  to  itself  in  the  modified  salesman 
problem,  the  size  of  the  solution  set  is  smaller  than  n!,  the  usual 
size  for  assignment  problems.   The  table  below  illustrates  solution- 
set  sizes  for  several  selected  problems: 
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Number 
of  cities 
(n) 


Number  of 
Assignment 
Solutions 
(n!) 

2 

6 

24 

120 


Number  of 
Modified  Traveling- 
Salesman  Solutions 
n!  >  (?)£■  (n-1)! 

1 

2 

9 

44 


Number  of 
Solutions  that 
are  Tours 
(n-l)i 

1 

2 

6 

24 


A  generalized  expression  for  the  third  column  above  was  not  derived, 
but  it  is  clear  that  the  fraction  of  the  time  a  modified  salesman 
solution  is  also  a  tour  is  greater  than  (n-l)!/n!  =  1/n. 

An  optimal  solution  to  the  modified  salesman  problem  provides  a 
lower  bound  on  the  minimum-distance  solution  to  the  original  problem. 
This  is  true,  because  the  modified  problem  differs  from  the  original 
only  in  that  it  has  fewer  constraints.   If  the  optimal  assignment 
solution  to  the  modified  problem  provides  a  tour,  then  the  original 
problem  is  solved;  if  not,  further  steps  must  be  taken. 

To  express  the  assignment  formulation  in  the  notation  of  Chapter  I, 
let  x. .  =  1  or  0  according  to  whether  or  not  the  salesman  is  to  travel 
from  city  i  to  city  j;  i,j  =  l,...,n.   Let  c..  be  the  distance  between 
city  i  and  j  for  all  i  and  j.   Set  p..  arbitrarily  large  for  all  i. 
Then, 


minimize 


T"   —  c    x 
1    j  ■  "ij  ij 


subject  to 


i    Xii  =  1   (J  =  1>-">n) 


and 


and 


-7—  x   =1   (i=l,.  .*  ,n) 
1     ij 


x.  ^  0  for  all  i,j 
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Alternately,  this  formulation  can  be  expressed  in  matrix  notation 

as  follows: 

T 
minimize  c  x  =  z 

subject  to  Ax  =  b 

and  x  s£  0 


where  the  definitions  of  c  ,x,b,  and  A  can  be  indicated  most  con- 
veniently by  example.   Consider  the  three-city  problem: 

T 
c    (cn,  c12,  c13,  c21,  c22,  c23,  c31,  c32,  c33), 


xll 

1 

x12 

1 

x13 

1 

X21 

1 

X 

X 

22 
x23 
x31 
x32 
x33 

> 

b 

1 
1 
1 
1 
1 

J 

1 

1 

1 

0 

0 

0 

( 

)     0 

0 

0 

0 

0 

1 

1 

1 

( 

)     0 

0 

0 

0 

0 

0 

0 

0 

] 

L     1 

1 

and  A  = 

1 

0 

0 

1 

0 

0 

] 

L     0 

0 

0 

1 

0 

0 

1 

0 

( 

)     1 

0 

0 

0 

1 

0 

0 

1 

( 

)     0 

1 

T         2        2  2  2 

In  general,  c   is  1  x  n  ,  x  is  n  x  1,  b  is  n  x  1,  and  A  is  2n  x  n  . 

Notice  that  since  the  k   column  of  A  corresponds  to  the  k    x,  .  of  x. 
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th  Nth  ,  th 

the  i   element  and  the  (n+j)    element  of  the  k   column  of  A  equal 

one  and  all  other  elements  in  that  column  equal  zero. 

The  assignment  problem  has  several  important  features  that  are 

listed  below  without  proof.   (see  references  3  and  5  for  proofs.) 

1.  The  assignment  problem  has  a  feasible  solution, 

2.  The  rank  of  the  matrix  A  is  2n-l,  thus  any  row  may  be  deleted 
as  redundant. 

3.  Any  set  of  2n-l  linearly  independent  column  vectors  from  A 
forms  a  basis  B. 

4.  The  inverse  of  any  basis  B  is  composed  entirely  of  zero,  one 
and  minus-one  elements. 

5.  Any  basic  feasible  solution  to  Ax  =  b  is  degenerate.   In 
fact,  there  are  exactly  n  basic  variables  equal  to  one  and  exactly 
n-1  equal  to  zero,  assuming  all  nonbasic  variables  are  arbitrarily 
set  equal  to  zero. 

C.   THE  SIMPLEX  SOLUTION  OF  THE  ASSIGNMENT  FORMULATION 

There  are  alternate  methods  for  solving  the  assignment  problem  but 
the  following  discussion  assumes  some  variant  of  the  simplex  method  is 
chosen.   The  final  simplex  tableau  contains  not  only  the  optimal  solu- 
tion to  the  assignment  problem,  but  other  interesting  information,  as 
well.   Before  proceeding,  a  few  more  notational  conventions  and 
definitions  are  stated  below: 

1.  Redefine  the  coefficient  matrix  A  such  that  one  row  (any  row) 
is  deleted. 

2.  Temporarily  redefine  the  vector  x  such  that  the  subscripts  on 

2 
its  elements  run  from  1  to  n  rather  than  1,1  to  n,n. 
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3.  Make  the  same  change  in  subscript  notation  for  the  elements 

of  c  . 

th  2 

4.  Denote  the  j    column  of  A  as  a.  for  all  j  =  1,...,  n  . 

5.  Let  B  denote  a  basis  in  A,  that  is,  any  set  of  2n-l  column 
vectors/ in  A. 

6.  Let  Xt,  denote  the  vector  of  x.'s  corresponding  to  the  a  '  s  in 

B  J  6         j 

B.   Call  the  elements  of  x™  the  values  of  the  basic  variables.   Note, 

since  xR  =  B  b,  x  is  simply  the  vector  obtained  by  summing  each  row 
B  B 

of  B.   Requiring  xR>  0,  then,  insures  its  elements  are  ones  and  zeros, 

exclusively. 

T 

7.  Let  c  denote  the  vector  of  c.'s  corresponding  to  the  elements 

B  J 

in  Xg. 

-1 

8.  Define  y.  =  B  a..   Note  that  the  y.'s  are  composed  entirely 

-*       J  J 

of  zero,  one,  and  minus-one  elements. 

T 

9.  Define  the  relative  cost  for  x.  as  c .  -  c  -c^y  . 

J     J    j   B  j 

The  simplex  procedure  involves  finding  a  basis  B  such  that  x  = 

B 

B  b  minimizes  z.   A  feasible  solution  is  recognized  when  the  elements 

of  xg  are  zeros  and  ones,  exclusively,  and  an  optimal  solution,  when 

>  2 

Ci  —   0  for  all  j  =  l,,,.,n  . 

Notice  that  c.  =  0  always  for  all  x.'s  in  xD  because  for  such 
J         y  J       B 

T 
x.  s,  c.  =  cu  y . .   If  c.  =  0  at  optimality  for  a  nonbasic  x.,  then  an 

alternate  optimum  can  be  obtained  by  bringing  that  x.  into  the  basis. 

In  this  case,  the  simplex  methodology  specifies  which  basic  variable 

is  to  be  replaced  in  order  to  maintain  feasibility. 

To  appreciate  the  impact  of  the  relative  costs  on  the  value  of  z, 

consider  the  following  relationship: 

z  =  z  +  ^~  c  .x. 
J€R  J  J 
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A 
where  z  is  the  new  value  of  z  and  R  is  the  index  set  of  all  nonbasic 

variables.   This  equation  shows  that  if  nonbasic  variables  are  set 
equal  to  one,  z  is  increased  by  the  sum  of  their  relative  costs. 

Of  course,  if  one  or  more  nonbasic  variables  are  set  equal  to  one, 
some  change  must  take  place  in  the  Xg  vector.   This  is  true  since  an 
assignment  solution  requires  exactly  n  variables  equal  to  one,  includ- 
ing the  nonbasic  variables.   In  order  to  measure  such  effects  on  Xg., 
consider  the  following  equality: 

A  "> 

x„  =  x  -  ^—    y  x 

B    B    j&R  'j  j 
A 

where  x  is  the  new  set  of  values  for  x„ .   Normally  x  =  0  for  all  i  in 

B  B  y      j  J 

R,  but  if  one  or  more  such  nonbasic  variables  are  set  equal  to  one,  the 
relation  above  changes  the  values  of  the  basic  variables.   Setting  one 

or  more  x.'s  not  in  x  equal  to  one  is  said  to  be  constraint  feasible 

J  B   H 

as  long  as  the  elements  of  x  are  zeros  and  ones,  exclusively. 

B 

D.   STATEMENT  OF  THE  ALGORITHM 

The  algorithm  under  study  begins  with  the  optimal  assignment 
solution  to  the  modified  traveling- salesman  problem,   If  the  optimal 
assignment  solution  provides  a  tour,  then  the  original  problem  is 
solved.   Otherwise,  the  algorithm  begins  to  introduce  combinations  of 
nonbasic  variables  at  the  one  level  until  a  tour  is  achieved. 
Complexities  arise  in  the  measures  taken  to  assure  that 
combinations  are  considered  in  ascending  order  of  their  relative  costs. 
If  these  steps  are  not  taken,  there  is  a  possibility  that  nonbasic 
variables  not  yet  tried  yield  lower-cost  tours. 

In  the  light  of  the  background  developed  in  the  previous  sections, 
the  algorithm  is  as  follows: 
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1.  Drop  the  constraint  that  a  traveling- salesman  solution  must 
provide  a  tour. 

2.  Solve  this  modified  problem  as  an  assignment  problem  using  a 
variant  of  the  simplex  method. 

3.  Check  the  optimal  assignment  solution  for  a  tour.   If  a  tour 
is  provided,  stop.   Otherwise,  continue. 

4.  Drop  simplex  tableau  vectors  corresponding  to  the  original  x. . 

variables.   Drop  tableau  vectors  corresponding  to  basic  variables. 

2  2 

There  are  n  -3n+l  vectors  remaining.   Let  m  =  n  -3n+l. 

5.  Rank  tableau  vectors  in  ascending  order  of  their  respective 
relative  costs. 

6.  Let  x-  =  x.,y.  =  v.,  and  c.  =  c.  where  i  =  1  if  c.  is  first 

J    i  j    i'      j    i  J 

in  the  ranking,  i  =  2  if  c-  is  second,  and  so  forth,  through  i  =  m  if 
c-  is  last  in  the  ranking.  (Retain  information  required  to  reconvert 
subscripts  as  needed  when  checking  for  tours.) 

7.  Let  i  =  0. 

8.  Increment  i  by  one. 

9.  If  x.  =  1  is  not  constraint  feasible  and  i  =  1,  return  to 
step  8.   If  X|  =  1  is  not  constraint  feasible  and  i>  1,  go  to  step  11. 
If  x^  =  1  is  constraint  feasible,  continue. 

10.  If  x.  =1  and  x_,  =  x  -y  provide  a  tour,  go  to  step  21.   Other- 

l  d  r>   i 

wise,  return  to  step  8  if  i  =  1;  continue  if  i  >  1. 

11.  Generate  all  combinations  of  x.  =1  with  x  =  1,  x  =1,..., 

i  12 

x.  i  =  1.   (Note  this  is  construed  to  mean  x.  =  1  appears  in  each 
combination  along  with  one  or  more  of  the  other  variables  indicated  set 
equal  to  one.   There  are  2    -  1  such  combinations.) 
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12.  Disregard  all  combinations  with  relative  costs  greater  than 
c.+1.   (Define  c^  ±   -     <**  . ) 

13.  Disregard  all  remaining  combinations  that  are  not  constraint 
feasible. 

14.  Check  remaining  combinations,  in  ascending  order  of  their 
relative  costs,  for  tours.   As  soon  as  a  combination  is  found  to  pro- 
vide a  tour,  go  to  step  21.   If  no  combination  yields  a  tour,  continue. 

15.  If  i  =  2  and/or  c.  ,  =  £•,  return  to  step  8.   Otherwise 
continue. 

16.  Let  k  =  i. 

17.  Decrement  k  by  one. 

18.  Generate  all  combinations  of  x.  =  1  with  x1  =  1,  X2  =  1,..., 

19.  Same  as  steps  13,  14,  and  15. 

20.  If  k  >  2,  return  to  step  17.   If  k  =  2,  return  to  step  8. 

21.  If  coming  from  step  19,  let  i  =  k.   Otherwise  do  not  change  i. 

22.  Compute  total  distance  for  the  tour  found.   Letting  K  be  the 

index  set  of  nonbasic  variables  included  in  the  solution,  the  total 

distance  for  this  tour  is  z  =  z  +~T„   c,  , 

k4K  k 

23.  If  i  ^  2,  go  to  step  32.   Otherwise,  continue. 

24.  If  all  combinations  of  nonbasic  variables  preceding  x.  have 
already  been  checked  against  the  current  value  of  z-z  as  an  upper 
bound,  go  to  step  32,   Otherwise,  continue. 

25.  Decrement  i  by  one. 

26.  Generate  all  combinations  of  x.  =  1  with  x..  =  1,  x  =  1,..., 


i  1    '   2 


Vl  =  l 


A 

z-z. 


27.   Disregard  all  combinations  with  relative  costs  greater  than 
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28.  Disregard  all  remaining  combinations  that  are  not  constraint 
feasible. 

29.  Check  remaining  combinations,  in  ascending  order  of  their 
relative  costs,  for  tours.   As  soon  as  a  combination  is  found  to  pro- 
vide a  tour,  go  to  step  30.   If  no  combination  yields  a  tour,  go  to 
step  31. 

30.  Update  the  value  of  z  for  the  latest  tour  found  and  continue. 

31.  If  i  y   2,  return  to  step  25.   Otherwise,  continue. 

32.  The  latest  tour  found  and  its  associated  z  is  an  optimal 
traveling-salesman  solution.   Stop. 

Since  eventually  all  combinations  of  nonbasic  variables  are  con- 
sidered, the  algorithm  must  eventually  find  a  tour.   Suppose  a  tour 
is  first  found  using  x.  =  1  in  some  combination  with  the  preceding 
nonbasic  variables.   This  tour's  relative  cost  is  less  than  or  equal 
to  c   -  because  only  combinations  meeting  this  criterion  were  considered, 
Thus,  there  is  no  need  to  try  combinations  beyond  x. .   There  is  a  need, 
however,  to  reconsider  combinations  prior  to  x. .   For  example,  x  =  1 
in  combination  with  x,  =  1  may  have  been  rejected  earlier  because  its 
relative  cost  is  greater  than  C-,  but  it  may  be  less  than  z-z.   Thus, 
proceeding  backward  through  the  nonbasic  variables  using  the  current 
value  of  z-z  as  an  upper  bound  guarantees  the  minimum  z:  will  be  found. 


21 


III.   THE  COMPUTER  PROGRAM 

The  algorithm  was  programmed  for  the  IBM  System/360  Model  67 
Operation  System  using  FORTRAN  IV  (H  level)  programming  language 
/Ref.  7/.   A  listing  of  the  complete  program  is  provided  in  Appendix 
A. 

Several  routines  used  in  the  program  were  not  developed  as  part 
of  this  study  but  were  modified  for  use  here  as  needed.   Subroutines 
ASSIGN,  LIP,  ELEM,  and  PIVOT  were  supplied  by  Professor  Greenberg  and 
are  used  to  solve  the  initial  assignment  problem.   Subroutine  RAND  is 
a  simple  routine  for  generating  city-to-city  distances,  uniformly 
distributed  from  0  to  100.   The  distance  from  each  city  to  itself  is 
arbitrarily  set  equal  to  999.   Subroutine  TIMEIT  is  used  to  measure 
algorithm  execution  times  in  seconds.   Part  of  TIMEIT  appears  in 
assembly  language. 

The  following  list  gives  a  brief  description  of  the  various 
routines  developed  especially  for  this  study. 

1.  The  Main  Program  follows  the  sequence  of  algorithm  steps  as 
stated  in  Chapter  II.   Provision  is  made  to  either  read  in  a  cost 
matrix  or  to  generate  one  randomly.   Other  provisions  are  made  to 
abort  if  something  unexpected  goes  wrong  and  to  begin  execution  on 
the  next  problem. 

2.  Subroutine  PRINTC  writes  out  the  cost  matrix.   Also,  alternate 
storage  space  is  provided  for  the  costs  because  they  get  overlayed 

in  the  assignment  routines  by  the  relative  costs. 

3.  Subroutine  PRINTA  writes  the  optimal  assignment  solution  In 
tableau  form.   If  the  problem  is  larger  than  seven  cities  they  y. 
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column  vectors  are  not  printed  because  there  are  too  many  to  fit  on  a 
page.   Subroutine  PRINTA  is  also  used  again  to  print  the  optimal  assign- 
ment solution  after  artificial  variables  remaining  in  the  basis  have 
been  replaced. 

4.  Subroutine  SOLCHK  checks  constraint  feasible  solutions  for 
tours.   If  a  tour  is  found,  it  is  printed. 

5.  Subroutine  RANKC  ranks  the  relative  costs  in  ascending  order  of 
magnitude.   Actually,  the  costs  are  not  moved  around  in  storage,  but 
second-order  subscripts  are  assigned  to  indicate  the  desired  ranking. 

6.  Subroutine  REPART  was  devised  to  replace  any  artificial 
variables  left  in  the  assignment  solution  and  to  select  a  redundant 
row  for  deletion  from  the  final  tableau. 

7.  Subroutine  LINCOM  is  the  most  important  subroutine  in  terms 
of  the  algorithm  at  hand.   This  subroutine  generates  the  combinations 
of  the  nonbasic  variables,  compares  the  relative  cost  of  each  combi- 
nation with  an  upper  bound,  and  checks  for  constraint  feasibility. 

A.   DEFINITIONS 

Three  common  blocks  (storage  areas)  are  used  in  the  program. 
Common  Block  ASGBLK  is  used  exclusively  with  the  four  assignment 
routines.   Common  Block  SHARE  is  used  to  pass  information  from  the 
assignment  routines  and  from  Subroutine  RAND  to  the  traveling- sales- 
man algorithm  routines.   Common  Block  TSBLK  is  used  to  pass  values 
among  the  algorithm  routines.   For  easy  reference,  this  section 
defines  the  variables  in  common  in  the  order  of  their  appearance. 
Variables  in  ASGBLK  are  not  listed  because  they  are  only  needed 
within  the  four  assignment  routines  which  are  not  part  of  this 
study. 
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1.   Common  Block  TSBLK 

a.  COST  is  the  sum  of  the  relative  costs  included  in  a 
particular  combination  of  nonbasic  variables,  each  taken  at  the  one 
level. 

b.  CSTLC(1600)  is  initially  used  as  alternate  storage  for 

the  city-to-city  costs.   After  artificial  variables  have  been  replaced, 
it  is  used  to  store  relative  costs  of  combinations  in  ascending  order 
of  magnitude. 

c.  NLC(IOOO)  is  used  to  store  a  combination  identifying 
number  corresponding  to  each  value  in  CSTLC(1600).   With  these  numbers, 
combinations  can  be  regenerated  as  needed. 

d.  IC(1600)  is  used  to  store  the  subscripts  of  the  variables 
associated  with  the  ranked,  relative  costs.   For  example,  the  first 
value  in  IC(1600)  is  the  subscript  of  the  nonbasic  variable  with  the 
smallest  relative  cost. 

e.  KC(40)  is  used  in  Subroutine  PRINTC  to  write  out  the  city- 
to-city  costs  as  integers. 

f.  KB(120)  is  used  to  store  the  new  values  of  the  basic 
variables  when  one  or  more  nonbasic  variables  are  set  equal  to  one. 
The  one-values  for  the  nonbasic  variables  are  augmented  to  this  vector 
starting  after  the  last  basic-variable  value,  (Also,  see  description 
of  IB(120)  below.) 

g.  IAX(40)  is  used  in  Subroutine  PRINTA  to  write  out  a  row 
of  the  final  assignment  tableau. 

h.   IR0W(40)  and  ICOL(40)  are  used  in  Subroutine  SOLCHK  to 
store  the  first  and  second  subscripts,  respectively,  of  variables  after 
these  have  been  reconverted  from  single- subscript  form. 
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i.   LC(40)  is  used  to  store  a  binary  number  up  to  40  digits  in 
length.   If  the  k   digit  is  a"  one,  then  the  k   ranked,  nonbasic 
variable  is  to  be  included  in  the  particular  combination  at  the  one 

level.   If  the  k   digit  is  zero-^  then  the  k   ranked,  nonbasic  variable 

40 
is  not  to  be  included.   LC(40)  can  be  used  to  represent  up  to  2    -  1 

different  combinations,  one  at  a  time. 

j.   KZ  is  the  truncated  value  of  Z  in  Subroutine  PRINTA.   In 
Subroutine  SOLCHK,  KZ  is  the  truncated  value  of  XZ  =  Z  +  COST.   If  the 
city-to-city  costs  are  integers,  KZ  =  XZ  always. 

k.   NC  =  N-3MM+1  is  the  number  of  nonbasic  variables,  excluding 
xii  variaDl-es  and  artificials. 

1.   J  is  used  in  the  Main  Program  to  index  the  progression 
through  the  ranked,  nonbasic  variables. 

m.   ICJ  =  IC(J). 

n.   LI  -  LEG(1,ICJ)  and  L2  =  LEG(2,ICJ). 

o.   Jl  =  J-l. 

p.  CST  is  the  current  upper  bound  against  which  combinations 
are  compared.  Any  combination  of  nonbasic  variables  at  the  one  level 
with  a  relative  cost  greater  than  CST  is  disregarded. 

q.   XZ  =  Z+COST.   Subroutine  SOLCHK  truncates  this  value  when 
writing  it  out.   (See  definition  (j)  above.) 
2.   Common  Block  SHARE 

a.  C(1600)  is  used  initially  to  store  the  city-to-city  costs. 
After  the  assignment  solution,  the  relative  costs  are  stored  here. 
(Also,  see  description  of  CSTLC  (1600)  above.) 

b.  B(80,80)  is  the  inverse  of  the  optimal  assignment  basis. 
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c,  H(80)  is  used  to  store  the  values  of  the  basic  variables 
in  the  optimal  assignment  solution.   (Also,  see  description  of  KB(120) 
above . ) 

d.  Z  is  the  optimal  assignment  total  cost. 

e„  LEG(2,1600)  is  a  matrix  of  values  used  to  convert  single- 
subscripted  variables  back  to  their  original  double-subscripted  form. 
For  example,  x,  =  x. .  if  i  =  LEG(l,k)  and  j  =  LEG(2,k)  -  MM  where  MM 

is  the  number  of  cities.   Furthermore,  LEG  provides  a  convenient  device 

-1  th 

for  obtaining  the  product  B   a,  =  y,  .   It  is  known  that  the  i   and 

th 
(j+MM)    elements  of  a,  are  ones,  if  i  and  j  are  defined  as  above,  and 

all  other  elements  of  a   are  zeros.   Thus,  rather  than  multiplying 

B   times  a  ,  simply  add  the  i   and  (j+MM)    elements  of  each  row  of 

-1 

B   to  produce  the  vector  y,  „ 

f.  IB(120)  is  used  to  store  the  subscripts  of  the  optimal 
assignment  basic  variables.   Negative  subscripts  indicate  artificial 
variables.   Starting  just  after  the  last  basic-variable  subscript, 
nonbasic-variable  subscripts  are  appended  to  this  vector  for  each 
nonbasic  variable  set  equal  to  one  in  a  particular  combination.   (Also, 
see  description  of  KB(120)  above.) 

g.  MM  is  the  size  of  the  problem,  that  is,  the  number  of 
cities. 

h.   M  is  the  number  of  constraining  equations  in  the  assignment 
formulation.   Initially,  M  =  2 (MM) ,  but  in  Subroutine  REPART  it  is 

reduced  by  one. 

2 
i.   N  =  (MM)   is  the  total  number  of  variables,  excluding 

artificials. 
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B.   DESCRIPTION  OF  SELECTED  SUBROUTINES 

The  Main  Program  was  written  to  conform  closely  to  the  algorithm 
as  presented  in  Chapter  II  and  can  be  followed  readily  from  the  listing 
in  Appendix  A.   Further  amplification  of  Subroutines  LINCOM,  SOLCHK, 
and  REPART  is  considered  necessary,  however,  in  order  to  clarify  the 
workings  of  the  overall  program.   Subroutine  REPART  was  devised  only 
because- the  assignment  routines  used  leave  artificials  in  the  basis 
and  do  not  drop  a  redundant  equation.   Subroutines  SOLCHK  and  LINCOM 
represent  techniques  created  to  fulfill  algorithm  requirements  and,  as 
such,  are  the  most  important  products  of  this  study. 

1.   Subroutine  REPART 

Subroutine  REPART  scans  the  subscripts  of  the  basic  variables 
returned  from  the  assignment  routines.   As  soon  as  a  negative  sub- 
script is  located,  its  corresponding  artificial  variable  is  replaced 
in  the  basis  by. ..the  first,  ranked,  nonbasic  variable  that  has  a  nonzero 
element  in  the  corresponding  row  of  its  y.  vector.   Using  this  nonzero 
element  as  a  pivot  element,  standard  row  transformations  are  performed. 
The  result  is  that  the  nonbasic  variable  is  brought  into  the  basis  and 
the  artificial  is  dropped.   The  optimal  value  of  the  objective  function 
does  not  change  because  the  artificial  variable  was  necessarily  at  the 
zero  level  implying  its  replacement  must  also  be  at  the  zero  level. 
After  replacing  an  artificial,  relative  costs  for  the  remaining  non- 
basic  variables  are  recomputed  and  these  variables  are  then  ranked  in 
ascending  order  of  their  new  relative  costs.   The  above  process  is 
repeated  until  all  artificials  have  been  removed  from  the  basis. 

It  is  known  that  one  of  the  assignment- tableau  rows  is  redundant, 
If,  while  trying  to  replace  an  artificial,  a  nonzero  pivot  element 
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cannot  be  found,  that  row  is  earmarked  for  deletion.   Otherwise,  after 
all  artificials  have  been  replaced,  each  row  is  scanned  for  all  zero 
elements.   If  such  a  row  is  found,  it  is  dropped.   Otherwise,  the  last 
tableau  row  is  arbitrarily  selected  for  deletion.   If  other  than  the 
last  row  is  dropped,  remaining  rows  are  shifted  upward  in  the  tableau 
to  fill  the  gap. 

Normally,  control  is  returned  to  the  statement  following  the 
calling  statement  in  the  Main  Program.   If  for  any  reason,  however, 
more  than  one  row  has  gone  to  all  zeros,  control  is  returned  to  state- 
ment 10  in  Main  and  execution  on  the  next  problem  is  started. 
2.   Subroutine  SOLCHK 

Subroutine  SOLCHK  examines  an  assignment  solution  to  see  if  it 
provides  a  tour,   Its  approach  is  to  count  the  number  of  steps  required 
to  return  to  the  city  of  origin.   If  this  number  is  less  than  the  total 
number  of  cities  MM,  then  subtours  exists  and  control  is  returned  to 
the  statement  following  the  calling  statement  in  Main,   If  the  number 
of  steps  required  is  exactly  equal  to  MM,  the  solution  is  a  tour.   In 
the  latter  event,  the  solution  is  printed  with  its  associated  total 
distance  and  control  is  returned  to  the  statement  in  Main  indicated  in 
the  parenthesis  of  the  calling  statement. 

To  illustrate  the  bookkeeping  involved,  consider  the  two  fol- 
lowing four-city  assignment  solutions: 


Steps 


x8  =  x2,4  =  l  l 

x3  =  xl,3  =  l 

x14  =  x4,2  =  l  2 

x9  =  x3,l  =  l 
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Steps 
X10  =  X3)2  =  l  l 

X13  =  X4,l  =  l  3 

X8  =  x2,4  =  l  2 

x2  =  x1>2  =  1  4 

Both  are  assignment  solutions,  but  only  the  second  provides  a  tour. 

Notice  that  the  routine  will  keep  scanning  the  KB  vector  until 
the  (MM)    one-value  is  found  and  then  scanning  is  terminated.   This 
approach  is  justified  because  only  constraint  feasible  solutions  are 
sent  to  SOLCHK  and  these  solutions  necessarily  have  exactly  MM 
variables  equal  to  one,  including  any  nonbasic  variables  set  equal 
to  one. 

3.   Subroutine  LINCOM 

The  essence  of  the  algorithm  lies  in  altering  the  optimal 
assignment  solution  by  setting  selected  nonbasic  variables  equal  to 
one  until  a  tour  is  achieved.   Subroutine  LINCOM  was  devised  to  generate 
all  the  required  combinations  of  nonbasic  variables.   An  attempt  was 
made  to  program  LINCOM  so  that  the  number  of  combinations  necessary  to 
consider  be  held  to  a  minimum.   The  objective  was  to  generate  them  in 
ascending  order  of  total  relative  costs  so  that  as  soon  as  a  combina- 
tion's cost  exceeded  the  current  upper  bound,  no  further  combinations 
need  be  considered.   This  objective  was  only  partially  achieved. 

The  most  difficult  problem  lay  in  the  fact  that  there  is  no 
way  to  completely  predetermine  the  ranking  of  all  possible  combinations. 
Consider  the  case  where  Ci  ^  ^2  —  •■••    cq—  ^fi  anc^  x_  =  1  is  found  not 
constraint  feasible.   Before  testing  xfi  =  1,  all  combinations  of 
Xr  =  1  with  the  preceding  nonbasic  variables  must  be  tried  because  some 
may  have  total  relative  costs  less  than  or  equal  to  c,. 
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One  way  to  test  the  combinations  would  be  to  try  x,  =  1,  x-  =  1. 
and  x,  .=.  1  singly  with  Xr  =  1,  then  in  pairs,  then  triplets,  and 
finally  all  four.   A  better  approach,  however,  is  to  set  up  the  four- 
digit  binary  number  0000  and  then  simply  increment  it  by  one,  fifteen 
times.   Number  the  digit  positions  from  one  through  four,  from  right  to 
left.   Then,  if  a  one  appears  in  the  j   position  where  j  =  1,...,4, 
let  x.  =  1  appear  in  the  combination.   The  possible  combinations  can 
now  be  represented  as  shown  below: 

Identifying 

Number  Combination 

1  0001 

2  0010 

3  0011 

4  0100 

5  0101 

6  0110 

7  0111 

8  1000 

9  1001 

10  1010 

11  1011 

12  1100 
'  13  1101 

14  1110 

15  1111 

Combination  number  9,  for  example,  represents  x,.  =  1  with  x,  =  1  and 
xi  =  1  for  a  total  relative  cost  of  c^  +  c,  +  c^. 

Notice  that  the  following  generalizations  can  be  made  about 
the  respective  total  relative  costs  of  the  combinations  above: 

#1  ^  #2  <=.     #3, 
#4  ^   #5  <:  #6  ^  #7, 
#8  ^  #9  —  #10  ^  #11, 
and  #12  ^  #13  —  #14  ^  #15. 
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Furthermore,  #3  ^  #5  ^  #9,  but  #3 ' s  relationships  with  #4  and  #8 

cannot  be  stated,  in  general.   Also,  #7  ~  #11  —  #13,  but  #7 '  s 

relationships  with  #8,  #9,  #10,  and  #12  are  unknown.   Finally,  #11^: 

#13,  but  there  is  no  way  to  tell  in  advance  which  is  smaller,  #11  or 

#12. 

While  there  are  always  many  organizational  difficulties 

associated  with  writing  and  debugging  a  large  program,  the  situation 

discussed  above  presented  one  of  the  most  challenging  theoretical 

problems  to  resolve.   The  scheme  finally  chosen  is  stated  below: 

If  a  combination's  relative  cost  is  greater  than  the  current  upper 
bound  and  the  combination's  identifying  number  is  equal  to 

i 

1)  2  ,  where  i  is  any  nonnegative  integer,  terminate  genera- 
tion of  further  combinations  (because  all  their  relative 
costs  will  exceed  the  current  upper  bound),  or 

2)  2  +  2  ,  where  i  and  j  are  any  nonnegative  integers  such 
that  j  <  i,  skip  to  combination  number  21+1  and  continue 
generating,  or 

3)  otherwise,  continue  generating. 

In  most  cases,  the  above  procedure  greatly  reduces  the  number  of 

combinations  necessary  to  generate.   It  could  be  expanded  by  saying 

i  j  k 
that  if  the  combination  number  is  equal  to  2+2+2  ,  where  i,  j,  and 

k  are  nonnegative  integers  such  that  k<j<i,  skip  to  number  2  +2    , 

and  so  forth.   Except  for  very  large  problems,  however,  such  added 

checks  would  probably  consume  more  computer  time  than  they  would  save. 

Another  facet  of  the  problem  of  trying  to  generalize  the 

relationships  among  relative  costs  is  presented  when  a  combination  is 

found  constraint  feasible.   To  illustrate  the  situation,  suppose 

combination  number  11  is  found  to  provide  a  relative  cost  less  than  the 

current  upper  bound  and  is  also  found  constraint  feasible.   If  #11  is 

sent  to  SOLCHK  at  this  point  and  found  to  provide  a  tour,  there  can  be 
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no  certainty  that  it  is  a  least-cost  tour  among  the  current  set  of  all 
possible  combinations.   This  is  because  #12,  as  yet  unchecked,  may 
provide  an  even  smaller  relative  cost,  may  be  constraint  feasible,  and 
may  also  provide  a  tour. 

The  technique  devised  to  resolve  the  above  difficulty  was  to 
store  the  identifying  number  and  relative  cost  of  each  combination 
found  constraint  feasible.   With  each  addition  to  storage,  the  numbers 
and  relative  costs  are  reranked  in  ascending  order  of  relative  costs. 
Finally,  when  all  necessary  combinations  have  been  generated  and 
checked,  the  constraint  feasible  solutions  are  sent  to  SOLCHK  in 
ascending  order  of  relative  costs,  thereby  insuring  the  first  tour 
found  has  a  cost  less  than  or  equal  to  any  tour  the  remaining  solutions 
in  storage  might  provide.   As  soon  as  a  tour  is  found,  control  is 
returned  to  the  first  address  indicated  in  the  parenthesis  of  the 
calling  statement  in  Main,   If  none  of  the  solutions  in  storage  provides 
a  tour,  control  is  returned  to  the  statement  in  Main  following  the 
calling  statement. 

Finally,  a  word  about  the  regeneration  of  combinations  required 
by  the  above  procedure.   The  decimal  identifying  number  in  storage  is 
converted  to  binary  through  the  standard  technique  of  successive 
divisions  by  2.   Remainders  become  the  binary  digits  from  right  to  left. 
The  binary  number,  then,  is  the  original  representation  of  the  combina- 
tion represented  temporarily  in  storage  by  the  decimal  identifying 
number.   Using  the  binary  number,  the  combination  of  interest  is  re- 
constructed and  sent  to  SOLCHK  along  with  its  associated  relative 
cost. 
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C.   LIMITATIONS  OF  THE  PROGRAM 

The  storage  areas  of  the  program  are  set  up  to  deal  with  problems 
smaller  than  41  cities,  but  these  areas  can  be  expanded  easily  if 
necessary.   The  four  subroutines  used  to  provide  an  optimal  assignment 
solution  are  programmed  to  terminate  execution  in  the  unlikely  event 
that  more  than  299  iterations  are  required.   The  only  other  significant 
limitations  of  the  program  are  found  in  Subroutine  LINCOM,  but  one  of 
these  presents  a  serious  restriction. 

In  Subroutine  LINCOM,  up  to  1000  constraint  feasible  solutions  can 
be  stored  at  one  time,  but  this  does  not  appear  to  present  an  important 
limitation.   In  any  case,  storage  space  for  CSTLC  and  NLC  in  Common 
Block  TSBLK  can  be  increased  considerably  if  necessary.   The  most 
serious  restriction  and  the  one  hardest  to  correct  is  the  fact  that 
Subroutine  LINCOM  can  generate  all  the  possible  combinations  for  no 

more  than  32  nonbasic  variables.   This  means  that  if  a  tour  is  not 

rd 
found  upon  setting  the  33   ranked,  nonbasic  variable  equal  to  one, 

then  execution  on  this  problem  must  terminate.   The  reason  for  this 

situation  is  that  each  combination  must  have  a  unique  identifying  number 

associated  with  it  in  LINCOM  and  the  largest  integer  number  that  can 

31  32 

be  expressed  on  the  IBM  360  is  2   -1.   There  are  2   -1  possible 

rd 
combinations  of  the  33   variable  with  the  32  preceding  it,  hence, 

the  problem.   As  an  example  of  the  gravity  of  the  restriction,  note 
that  a  forty-city  problem  presents  1481  nonbasic  variables  for 
consideration. 

The  perplexity  above  is  not  insurmountable,  but  it  was  not  re- 
solved in  this  study.   The  possibility  of  representing  the  identifying 
number  in  floating-point  mode  is  precluded  because  a  form  of  the 
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number  is  used  as  a  DO  Loop  parameter.   (See  the  variable  KT  in  LINCOM) 

It  might  be  possible  to  represent  the  number  as  an  element  of  a  vector. 

31 
For  example,  each  of  the  first  2   -1  numbers  might  be  prefixed  with  a 

31 
1,  the  second  set  of  2   -1  numbers  with  a  2,  and  so  forth.   The  DO 

Loop  using  the  number  as  a  parameter  could,  in  like  manner,  be 

expanded  into  a  double  DO  Loop. 
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IV.   COMPUTATIONAL  EXPERIENCE 

Unfortunately,  time  constraints  on  this  study  did  not  permit 
extensive  testing  of  the  computer  program  in  its  final  form.   Enough 
was  accomplished,  however,  to  indicate  several  important  conclusions. 

Appendix  B  presents  a  partial  facsimile  of  the  program  output  for 
the  same  ten-city  problem  that  is  discussed  in  References  8  and  9. 
Eastman's  algorithm  required  the  solution  of  eleven  assignment  prob- 
lems before  the  optimal  tour  was  obtained.   Little's  method  required 
31  branches  to  arrive  at  the  desired  solution.   Appendix  B  indicates 
the  programmed  algorithm  under  study  took  only  1.50  seconds  after  the 
initial  assignment  solution  to  obtain  the  optimal  tour.   The  first 
tour  was  found  after  trying  combinations  of  the  first  23  ranked,  non- 
basic  variables  and  this  tour  turned  out  to  be  the  optimal  one. 

City-to-city  costs  were  generated  by  Subroutine  RAND  for  17 
twenty-city  problems.   Fourteen  of  these  problems  were  solved  by  the 
program  in  a  mean "time  of  25  seconds  and  after  trying  an  average  of  15 
nonbasic  variables  for  each  problem.   Execution  was  terminated  on  the 
other  three  problems  because  a  tour  had  not  been  found  within  the  first 
33  nonbasic  variables. 

Ten  randomly  generated  thirty-city  problems  were  attempted  and 
nine  were  solved  in  an  average  time  of  56  seconds.   The  mean  number 
of  nonbasic  variables  required  for  testing  in  these  nine  problems  was 
16. 

Finally,  Subroutine  RAND  was  used  to  create  three  forty-city 
problems.   The  first  ran  .07  seconds  after  the  assignment  solution  to 
discover  this  was  also  an  optimal  tour.   The  second  problem  ran  549 
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rd 
seconds  before  terminating  after  trial  of  the  33   nonbasic  variable 

failed  to  provide  a  tour.   The  third  problem  was  solved  in  187 

seconds  after  trying  only  four  nonbasic  variables. 

In  view  of  the  computational  experience  cited  above,  one  can 

conclude  that  the  algorithm  works  and  that  the  program  works  with  one 

exception.   Of  the  31  problems  tried,  five  could  not  be  solved  because 

tours  had  not  been  found  upon  setting  the  33   nonbasic  variables 

equal  to  one.   If  the  program  is  to  be  of  general  use,  this  somewhat 

arbitrary  restriction  in  Subroutine  LINCOM  must  be  removed.   The 

limitation  could  be  withdrawn  if  a  way  is  found  to  obviate  the 

identifying  numbers  presently  associated  with  each  possible  combina- 

31 
tion.   Since  2   -1  is  the  largest  integer  that  can  be  represented  on 

the  IBM  System/360,  none  beyond  the  32   ranked,  nonbasic  variable  can 

be  tried  in  combination  with  its  preceding  variables. 

Several  steps  can  be  taken  to  increase  the  execution  speed  of  the 
program.   All  write  statements,  except  those  in  Subroutine  SOLCHK  and 
the  various  termination  messages  throughout  the  program,  can  be  replaced 
with  CONTINUE  statements.   In  subroutine  REP ART  there  is  no  particular 
reason  why  relative  costs  are  recomputed  and  reranked  after  each 
artificial  replacement.   These  two  chores  could  be  postponed  until  all 
artificials  have  been  replaced.   Since  a  replacement  variable  enters 
the  basis  at  the  zero  level,  it  is  not  necessary  that  its  relative 
cost  be  the  least  possible. 

In  conclusion,  a  possible  improvement  to  the  algorithm,  itself,  is 
suggested.   The  algorithm  is  very  efficient  if  a  tour  is  found  within 
the  first  few  nonbasic  variables  tried,  otherwise  much  time  is  consumed 
rechecking  combinations  already  considered  against  a  new  upper  bound, 
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and  perhaps  again  against  an  even  newer  upper  bound,  and  so  forth. 
This  phenomenon  can  be  observed  in  the  example  output  displayed  in 
Appendix  B.   The  algorithm  might  be  improved  as  follows:   1)  obtain 
a  relatively  low-cost  tour  by  linking  together,  on  some  reasonable 
fashion,  the  subtours  appearing  in  the  optimal  assignment  solution, 
2)  compute  z  for  this  tour,  3)  using  z-z  as  an  upper  bound,  begin 
with  the  last,  ranked,  nonbasic  variable  whose  relative  cost  is  less 
than  z-z  and  proceed  backward  just  as  in  the  original  algorithm  after 
the  first  tour  is  found.   In  this  manner,  it  is  not  necessary  to 
proceed  forward  through  the  ranked,  nonbasic  variables  looking  for 
that  first  tour,  and  it  is  necessary  to  work  backwards  through  the 
variables  only  once. 
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APPENDIX  A    THE  PROGRAM  LISTING 


C   PLACE  GREEN  JOB  CARD  HERE. 

//   EXEC   FORTCALG, REGION. G0=160K, TIME. GO=10 

//FORT.SYSPRINT  DO  SPACE=( CYL, ( 4, 1) ) 

//FORT.SYSIN  DD  * 

C   A  TRAVELING-SALESMAN  ALGORITHM 

COMMON  /TSBLK/  CO  ST, CSTLC ( 1600 ) , NLC ( 1 000 ) , 

1  IC(  1600),KC(40),KB(  120)  , I  AX (40) , 

2  IR0W(40), ICOL(40) ,LC(40), 

3  KZ,NC, J,ICJ,Ll,L2, J1,CST,XZ 
COMMON  /SHARE/  C(1600),B(  80,  80), H(  80), Z, 

1  LEG(2, 1600) ,16(120), MM, M,N 

C   '4M  IS  THE  NUMBER  OF  CITIES. 

C   NRAND=1  IF  COSTS  ARE  TO  BE  RANDOMLY  GEMERATED. 
10  RFAD(5,5CC,END=90)  MM,NRAND 
500  FORMAT (21 2) 

M=2*MM 

N=MM*MM 

IF(NRAND.EQ.l)  CALL  RANDU20) 

READ(5,51C)  (C(J),J=1,N) 
510  F0RMAT(16F5.0) 
20  CALL  PRINTC 
C   OBTAIN  OPTIMAL  ASSIGNMENT  SOLUTION. 

CALL  ASSIGNU10) 

CALL  PRINTA 
C   START  CLOCK. 

CALL  TIMEIT(0,TIME) 
C   SEE  IF  OPTIMAL  ASSIGNMENT  SOLUTION  IS  A  TRAVELING- 
C   SALESMAN  SOLUTION.   IF  YES  GO  TO  80,  IF  NO  CONTINUE. 

COST=C. 

CALL  SOLCHKU80) 
C   RANK  NONBASIC  VECTORS  IN  FINAL  TABLEAU  IN  ASCENDING 
C   ORDER  OF  THEIR  ASSOCIATED  RELATIVE  COSTS. 

CALL  RANKC 
C   REPLACE  ARTIFICIALS  REMAINING  IN  OPTIMAL  ASSIGNMENT 
C   SOLUTION  AND  DROP  A  REDUNDANT  EQUATION  FROM  THE  TABLEAU. 
C   NOTE  THAT  M  IS  DECREMENTED  BY  1  IN  SUBROUTINE  REPART. 

CALL  REPART(SIO) 

CALL  PRINTA 

C(  IC(NC-H)  )  =  1000000. 

CSTHLD=C. 
C   NC=MM*MM-3*MM+1  IS  THE  TOTAL  NUMBER  OF  NONBASIC 
C   VARIABLES  OF  INTEREST. 
C   PROCEED  FORWARD  THROUGH  THE  RANKED,  NONBASIC  VARIABLES. 

DO  5C  L=1,NC 

J  =  L 

ICJ=IC(J) 

KCST=C(ICJ) 

WRITE(6,61C)  ICJ,KOST 
610  F0RMAT(//1CX, 'TRY  X(«,I2,»)  =  1    REL  COST  =  »,I3) 

L1=LFG(1,ICJ) 

L2=LEG(2,ICJ) 

IB(M+1  )  =  ICJ 

KP(M+1)=1 

COST  =  C(ICJ  ) 

DO  3C  1  =  1,  M 

XHOLD=H( I )-B(  I,L1)-B(  I,L2) 

KB (I )=XHOLD+SIGN( .5,XH0LD) 
C   IF  X(ICJ)=1  IS  NOT  CONSTRAINT  FEASIBLE,  GO  TO  40. 

IF(KB( I J.ME.O. AND.KB(  I  ).NE.l)  GO  TO  40 
30  CONTINUE 

CALL  S0LCHK(&58) 
40  IF(J.EQ.l)  GO  TO  50 
C   TRY  ALL  POSSIBLE  COMBINATIONS  OF  X(ICJ)=1  WITH  THE 
C   PRECEDING  NONBASIC  VARIABLES.   DISREGARD  ANY  COMBINATION 
C   WITH  RELATIVE  COST  GREATER  THAN  CST=C ( I C ( J+l )  )  . 

J1=J-1 

WPITE(6,620)  Jl 
620  F0RMATC10X, »TRY  THE  (2**,,I2,,>  -  1  LN  COM') 
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KJ1 


ICJ 
ICJ 


630 
IF 


CST=C(  IC ( J  +  l 
IF    A    COMBINATIO 
OTHERWISE    CONTI 
CALL    LINCOM( 
IF    CST    IS    EQUAL 
IF(CST.EQ.C( 
IF(J.EQ.2)    G 
PROCEEDING    BACK 
VARIABLES,    RECH 
CSTHLD=CST 
KJ1=J-1 
DO    45    K=2 
J=KJl-K+2 
J1=J-1 
ICJ=IC(J) 
Li-LEGilt 
L2=LEG(2, 
IB(M+1)=ICJ 
KB(M+1)=1 
WRITE<6,630) 
FCRMAT(//1^X 
WRITE(6,62C) 
A    COMBINATIO 
OTHERWISE    CONTI 
CALL    LINCOM< 
45    CONTINUE 
50    CONTINUE 
THE    PROGRAM    SHO 
THERE    ARE     (MM-1 
CALL    TIMEITI 
WRITE(6,64C) 
640    FORMAT(//10X 
1 
WPITE<6,65C) 
F0RMAT(//1CX 
GO    TO    IC 
CLC    COMBINAT 
OLD    CST    .LE. 
FOUND,    GO    TO    80 
58    IFCCST.LE.CS 
60    IF(J.LE.2)    G 
PROCEEDING    BACK 
VARIABLES,    RECH 
WHICH    IS    THE    RE 
KJ1=J-1 
DO    70    K=2 
J=KJl-K+2 
J1=J-1 
ICJ=IC( J) 
L1=LEG(1, 
L2=LEG(2, 
I6(M*1)=ICJ 
KB(M+1)=1 
WRITE(6,63C) 
WRITE(6,620) 
WHETHER    OR    NOT 
A    TOUR,     CONTINU 
CALL    LINCOM< 
70    CONTINUE 
STOP    CLOCK. 
80    CALL    TIMEIT( 
WRITE(6,66C) 
F0RMAT(//1CX 
WRITE(6,65C) 
GO    TO    IC 
STOP 
END 


)  ) 

N  PROVIDES  A  TOUR  GO  TO  58, 
NUE. 
&58,£ 
TO  T 
ICJ)  ) 
0  TO 
WARD 
ECK  0 


10) 

HE  PREVIOUS  CST  GO  TO  50. 

GO  TO  50 
50 

THROUGH  THE  RANKED,  NONBASIC 
LD  COMBINATIONS  AGAINST  THE  NEW  CST. 


650 

IF 

AN 


ICJ 

, 'RETURN  TO  X( •,  12,  ■  )  AND'  ) 

Jl 
N  PROVIDES  A  TOUR  GO  TO  60, 
NUE. 
£60, £101 


ULD  NEVER  REACH  THIS  POINT  BECAUSE 

):  POSSIBLE  TOURS. 

-1,TIME) 

,»A  TRAVELING  SALESMAN  SOLUTION1, 

•  DOES  NOT  EXIST* ) 

TIME 
.•ALGORITHM  TIME  =«,F10.4,»  SECONDS') 

IONS  HAVE  ALREADY  BEEN  CHECKED  AGAINST 
THE  RELATIVE  COST  OF  THE  TOUR  JUST 

THLD)  GO  TO  80 

0  TO  80 

WARD  THROUGH  THE  RANKED,  NONBASIC 

ECK  OLD  COMBINATIONS  AGAINST  THE  NEW  CST 

LATIVE  COST  OF  THE  LATEST  TOUR  FOUND. 


KJ1 


ICJ 
ICJ 
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90 


ICJ 
Jl 

A  COMBINATION  IS  FOUND  TO  PROVIDE 

E. 

£70, £10) 


-1,TIME) 

KZ 
,'**  Z  =,,I5,»  IS  OPTIMAL  **•) 

TIME 
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SUBROUTINE  PRI NTC 

COMMON  /TSBLK/  COST, CSTLC ( 1600 ) ,NLC ( 1000 ) , 

1  IC(1600),KC(40),KB( 120 ) . I  AX (40) , 

2  IR0W(40) ,ICOL(40) ,LC(40) , 

3  KZ,NC, J, ICJ,L1,L2, J1,CST,XZ 
COMMON  /SHARE/  C(1600),B(  80,  80), H(  80), Z, 

1  LEG(2,1600) ,IB(120) ,MM,M,N 

WRITE(6,6C0) 
6^0  FORMATC  1'  ,11X,'C0ST  MATRIX*) 

WRITE(6,610)  (I, 1=1, MM) 
610  F0RMAT(//11X,40I3> 

WRITE(6,620) 
620  FORMAT!'  •  ) 

1=0 

DO  2C  K=1,N,MM 

L=K+MM-1 

J  =  0 

DO  10  LL=K,L 

J  =  J+1 

KC(J)=C(LL) 
10  CONTINUE 

1=1  +  1 

foRITE(6,63C)  I  ,  (  KC  (  J  )  ,  J=l  ,  MM  ) 
630  FORMAT(  8X , I  2 , IX, 40  I  3 ) 
20  CONTINUE 

DO  30  J=1,N 
30  CSTLC(  J)=C(J) 

RETURN 

END 

SUBROUTINE  PRINTA 

COMMON  /TSBLK/  COST , CSTLC ( 1600 ), NLC ( 1000 ) , 

1  IC(  160  0) ,KC(40) ,KB( 120) , I  AX (40) , 

2  IROW(40), IC0L(40),LC(40), 

3  KZ,NC, J,ICJ,L1,L2, J1,CST,XZ 
COMMON  /SHARE/  C(1600),B(  80,  8H).H(  80), Z, 

1  LEG(2, 1600) ,IB(120),MM,M,N 

WRITE(6,600) 
600  FCRMATP1' ,10X, 'OPTIMAL  ASSIGNMENT  SOLUTION*) 
IF(M.E0. (2*MM) )  GO  TO  10 
WRITE(6,610) 
610  FCRMAT(19X,' (MODIFIED) • ) 
GO  TO  60 
10  KZ=Z 

DO  20  1=1, M 
20  KB(I )=H(I )+.5 
NC=0 

MMl=MM-»-l 
DO  50  J=1,N 
DO  3C  1=1, M 

IF(J.EQ.IBU) )  GO  TO  50 
30  CONTINUE 

DO  40  K=1,N,MM1 
IF(J.EQ.K)  GO  TO  50 
40  CONTINUE 
NC=NC+1 
IC(NC)=J 
50  CONTINUE 

60  IFCMM.LE.7)  GO  TO  80 
WRITE(6,62C) 
620  F0RMAT(//1CX, 'BASIS   B'//) 
DO  70  1=1, M 

WPITE(6f63C)  IB(I),  KB(I) 
630  F0RMAT(1CX,I4,2X, 12) 
70  CONTINUE 

WRITE(6,64C)  KZ 
640  F0RMAT(/15X, 'Z=', 14) 


40 


80 
650 

660 

670 


90 

680 
100 

110 

690 


•COEFF 
( ICU) 
BASI  S 


RETURN 

WRITE(6,650) 

FORMAT (//22X 

WRITE(6,660) 

F0RMAT(/10X. 

WRITE(6,67C) 

FORMAT(»  •) 

DO  ICO  1  =  1, M 

DO  90  J=1,NC 

XXHOLD=B( I.LEGt 1, IC( 

IAX< J)=XXHOLD+SIGN( . 

CONTINUE 

WRITE(6,68C)     IB(I 

FORMATdCX,  14, 2X, 

CONTINUE 

DO    11C    J=ltNC 

KC(J)=C(  ICU)  ) 

WRITE(6,690) 

FORMAT(/10X. 

WRITE(6,64C) 

RETURN 

ENC 


VECTORS' ) 

J=1,NC> 

B       ',3713) 


) ))+B(I ,LEG(2,IC( J) > ) 

,XXHOLD) 
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KB(I),     (IAX(J)  ,J=1,NC> 
2X,37I3) 


(KC( 
•REL 
KZ 


J) ,J=1,NC) 
COST*, 2X, 3713) 


SUBROUTINE     SOLCHM*) 

COMMON    /TSBLK/    CO  ST, C STLC ( 1600 ) ,NLC ( 1000 ) , 

1  IC(160C),KC(40),KB(120),I AX (41), 

2  IROW(40> ,IC0L(40) ,LC(40)  , 

3  KZ,NC, J, ICJ,L1,L2, J1,CST.XZ 
COMMON  /SHARE/  C(1600),B(  80,  80), H(  80), Z, 

1  LEG( 2, 1600)  ,IB<120)  ,MM,M,N 

PIPTS  OF  BASIC  VARIABLES  EQUAL  TO  ONE 
NAL  DOUBLE-SUBSCRIPTED  FORM. 


SUBSCI 
ORIGI 


10 


30 


40 
50 


610 


620 
60 

630 


L=1,MM 


C   CONVERT 

C       BACK    TO 
1=0 
80   20 
1  =  1+1 
IF(KB(  I  ).E 
IPOW(L)=LE 
ICOL(L)=LE 
20    CONTINUE 

C       COUNT    THE    NUM 

C   TO  GO  FROM  OR 
KSTEP=1 
K=l 

DO  40  L=2, 
IFUROW(L) 
KSTEP=KSTE 
IFUCOL(L) 
K=L 

GO    TO    30 
CONTINUE 
IF(KSTEP.L 
CST=COST 
XZ=Z+COST 
KZ  =  XZ 

WRITE(6,61 
F0RMAT(//1 
DO  60  L=l, 
WPITE(6,62 
F0RMAT(21X 
CONTINUE 
WRITE(6,63 
FCRMAT(/28 
RETURN  1 
END 


10 


Q.O)  GO  TO 
G(l,  IB(  I  )) 
G(2,  IB( I  ))-MM 

BER  OF  STEPS  REQUIRED 
IGIN  BACK  TO  ORIGIN. 


MM 

.NE.  ICOL(K)  )  GO  TO  40 

P+l 

.EQ.IROWI 1) )  GO  TO  50 


T.MM)  RETURN 


0) 

CX,«A  TRAVELING  SALESMAN  SOLUTION  IS 

MM 

0)  IROW(L) t ICOL(L) 

,»X(  •,  12,  •,  ',  12, •  )  =  !•) 


—  •//) 


C)  KZ 


Z  =  », 15) 


41 


SUBROUTINE  RANKC 

COMMON  /TSBLK/  COST, CSTLC ( 1600 ) ,NLC( 100C ) , 

1  IC  (  1600  ) , KC( 40 ) , KB ( 120 ) , I  AX ( 40 ) , 

2  IR0W(40) ,ICOL(40) ,LC(40), 

3  KZ,NC1J,ICJ,L1,L2! J1,CST.XZ 
COMMON  /SHARE/  C(l600),B(  80,  80), H(  80), Z, 

1  LEG<2,1600) ,IB(120> ,MM,M,N 

NC1=NC-1 

DO  20  K=1,NC1 

CMIN=C(IC(K) ) 

LMIN=K 

K1=K+1 

DO  IC  L=K1,NC 

IF(C( IC(L) ) .GE.CMIN)  GO  TO  10 

CMIN=C(IC< L) ) 

LMIN=L 
10  CONTINUE 

ICHOLD=IC( K) 

IC(K)=IC(LMIN) 

IC(LMIN)=ICHOLD 
20  CONTINUE 

RETURN 

END 

SUBROUTINE  REPART(«) 

COMMON  /TSBLK/  COST, CSTLC ( 1600 ), NLC ( 1000 ) , 

1  ICC  1600)  ,KC(40),K3U20),IAX(40), 

2  IROW(4G),ICOL(40) ,LC(40) , 

3  KZ,NC,J,ICJ,LI,L2,J1,CST,XZ 
COMMON  /SHARE/  C(1600),B(  80,  80). H(  80), Z, 

1  LEG(2, 1600), IBl 120) ,MM,M,N 

IIHOLD=0 

DO  120  11=1, M 
C   IB  IS  VECTOR  OF  BASIC-VARIABLE  SUBSCRIPTS. 
C   NEGATIVE  VALUES  INDICATE  ARTIFICIALS. 

IF(IB( II ).GT.O)  GO  TO  120 

DO  10  J=1,NC 

ICJ=IC( J) 

L1=LEG(1,ICJ> 

L2=LEG(2, ICJ) 

XHLD=B(II,L1)+B( I I,L2) 

IPIV=XHLD+SIGN< .5,XHLD) 

IF(IPIV.NE.O)     GO    TO    30 
10    CONTINUE 
C       IF    ONE    ARTIFICIAL    CAN    NOT    8E    REPLACED    WITH    A 
C       PIVOT    OPERATION,     EARMARK     ITS    ROW    FOR    DELETION. 
C       IF    THERE    IS    NO    PIVOT    ELEMENT    FOR    MORE    THAN    ONE 
C       ARTIFICIAL,    TERMINATE    EXECUTION    ON    THIS    PROBLEM. 

IFUIHOLD.NE.01    GO    TO    20 

IIHOLD=II 

GO  TO  12C 
20  WRITE(6,6C0) 
600  F0PMAT(//1CX,'CAN  NOT  REPLACE  ARTIFICIALS*) 

RETURN  1 
C   REPLACE  ARTIFICIAL  WITH  A  PIVOT  OPERATION. 
30  ie<II)=ICJ 

DO    5C     1=1, M 

IFd.EQ.II  )    GO    TO    50 

XHLD=B(  I,L1UB<  I,L2) 

IYIJ=XHLD+SIGN(.5,XHLD) 

IF( IYIJ.EQ.O)    GO    TO    50 

IQ=IYIJ/IPIV 

DO    40    L=1,M 
40    B(I,L)=B(  I  ,L)-IQ*BUI  ,L> 
50    CONTINUE 

IF(IPIV.EQ.l)     GO    TO    70 

DO    6C    L=1,M 
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60  B( II,L)=(-1. )*B(II,L) 

70  NC=NC-1 

IFCJ.EQ.NC4l)  GO  TO  90 
SHIFT  NONBASIC  VECTORS  LEFTWARD  TO  FILL  GAP. 

DO  80  L=J,NC 
80  IC(L)=IC(L4l) 
RECOMPUTE  RELATIVE  COSTS. 
90  DO  11C  J=1,NC 

ICJ=IC( J) 

L1=LEG(1,ICJ) 

L2=LEG(2,ICJ) 

CHOLD=CSTLC( ICJ) 

DO  100  1  =  1,  M 

IF(IBU).LT.O)  CX=10CC000. 

IF(  IB( I) .GT.O  )  CX  =  CSTLC( IB( I ) ) 
100  CHOLD=CHOLD-CX*(B<  I  ,L 1 ) +  B ( I , L2 ) ) 
110  C(  ICJ)=CHOLD 

CALL  RANKC 
120  CONTINUE 
FIND  A  REDUNDANT  TABLEAU  ROW  TO  DROP. 

M1  =  M 

M=M-1 

IFUIHOLD.NE.O)  GO  TO  150 

DO  14C  I  1  =  1, Ml 

DO  13C  J=1,NC 

ICJ=IC(J) 

L1=LEG(1, ICJ) 

L2=LEG(2,ICJ) 

XHLD=B(II,L1 )+B( II, L2) 

IYIJ=XHLD+SIGN(.5,XHLD) 

IF( IYI J.NE.O)  GO  TO  140 
130  CONTINUE 

IIHOLD=II 

GO  TO  15C 
140  CONTINUE 

RETURN 
150  IF( IIH0LD.FQ.M1)  RETURN 
SHIFT  TABLEAU  ROWS  UPWARD  TO  FILL  GAP. 

DO  170  1=1 IHOLD,M 

IB(I  )=IB(I+1) 

H(  n-HU  +  1  ) 

DO  160  J=1,M1 
160  B( I,J)=B(I+1,J  ) 
17C  CONTINUE 

RETURN 

END 

SUBROUTINE  LINCOM(*,*) 

COMMON  /TSBLK/  CO  ST, CSTLC ( 1600 ) ,NLC ( 1000  )  , 

1  IC(1600) ,KC(40),KB( 120 ), I  AX (40) , 

2  IR0W(40) ,IC0L(40) ,LC(40) , 

3  KZ,NC, J, ICJ,L1,L2, J1,CST,XZ 
COMMON  /SHA^E/  C(1600),B(  80,  80). H(  80), Z, 

1  LEG(2, 1600) ,IB( 120) ,MM,M,N 

IFU.LE.32  )  GO  TO  20 
WPITE(6,6CC) 
600  F0RMAT(15X,'CAN  NOT  PROCEED  BECAUSE  LN  COM 
1        'WILL  EXCEED  2**31  -  1' ) 
10  CALL  TIMEIK-1  ,TIME) 
WPITE(6,61C) 
610  F0RMATI//15X,' PROGRAM  LIMITATIONS  FORCE  ', 
1        'ALGORITHM  HALT* ) 
WPITE(6,62C)  TIME 
620  F0RMAT(/15X, 'ALGORITHM  TIME  =',F10.4,'  SECONDS') 
RETURN  2 
INITIALIZE  COUNTERS. 
20  NGEN=0 
NCSTF=0 
NCNF=C 
KOST=CST 
INITIALIZE  COMBINATION  GENERATER. 
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DO  3C  1=1, Jl 
30  LC(I  )  =  0 
C   BEGIN  GENERATING. 

C  Jl  IS  EQUAL  TO  J-l  AND  IS  THE  TOTAL  NUMBER  OF 
C  NONBASIC  VARIABLES  TO  BE  TRIED  IN  COMBINATION 
C   WITH  THE  J  TH  NONBASIC  VARIABLE. 

C   KT  IS  THE  NUMBER  OF  BINARY  NUMBERS  WITH  EXACTLY 
C   I  Jl  DIGITS  WHERE  THE  LEFTMOST  DIGIT  IS  ALWAYS  A  ONE. 

DO  17C  U1  =  1,J1 

KT  =  2**UJ1-1> 

DO  16C  11=1, KT 

DO  40  1  =  1,  Ul 

IF(LC(  D.EQ.O)  GO  TO  50 
40  LC(I)=0 
50  LC(I )=1 
C   INCREMENT  THE  NUMBER  OF  COMBINATIONS  GENERATED  BY  ONE. 

NGEN=NGEN+1 
C   COMPUTE  LATEST  COMBINATION'S  RELATIVE  COST. 

COST=C(ICJ) 

DO  60  1  =  1,  IJ1 

IF(LC<  D.EQ.l)  COST=COST+C(IC(I  )) 
60  CONTINUE 

IF(COST.LE.CST)  GO  TO  110 
C   THE  FOLLOWING  SEQUENCE  OF  STATEMENTS  DOWN  TO 
C   STATEMENT  11C  ARE  VARIOUS  CHECKS  TO  SEE  IF 
C   ENTIRE  BLOCKS  OF  COMBINATIONS  CAN  BE  SKIPPED. 

IF(II.GT.l)  GO  TO  70 

IF(NCNF.GT  .0  >  GO  TO  180 
65  WRITE(6,64C>  NGEN , KOST, NCSTF ,NCNF 
640  F0RMAT(15X,»NR  LN  COM  NECESSARY  TO  GENERATE  =»,I10, 

1  /15X,»NR  WITH  REL  COST  .LE.»,I4,«  =,,6X,U0, 

2  /15X,«NR  FOUND  CONSTRAINT  FEASIBLE  =,,3X,I10) 
RETURN 

70  DC  80  1=2,  Ul 

IF(II.EQ.2**( I-2)+l)  GO  TO  90 
80  CONTINUE 

GC  TO  16C 
90  DO  100  1  =  1, Ul 
100  LC(I )=1 

GO  TO  17C 
C   INCREMENT  BY  ONE  THE  NUMBER  OF  COMBINATIONS 
C   FOUND  LESS  THAN  OR  EQUAL  TO  THE  UPPER  BOUND. 

110  NCSTF=NCSTF+1 
C   CHECK  THE  LATEST  COMBINATION  FOR  CONSTRAINT  FEASIBILITY. 
DO  130  K=1,M 

XH0LD=H(K)-B(K,L1)-B(K,L21 
KB(K)=XHOLD+SIGN( .5,XH0LD) 
DC  12C  1  =  1, Ul 
IF(LC(  H.EQ.O)  GO  TO  120 

ici=icm 

XHCLD=KB(K)-B(K,LEG(1, ICI ) )-B( K, LEG( 2 , IC I ) ) 
KB(K)=XHOLD+SIGN( .5,XH0LD) 
120  CONTINUE 

IF(KB(K) .NE.O.AND.KB(K).NE.l)  GO  TO  160 
130  CONTINUE 
C   INCREMENT  BY  ONE  THE  NUMBER  OF  COMBINATIONS 
C   FOUND  CONSTRAINT  FEASIBLE. 
NCNF=NCNF+1 

IF  (NCNF.LE.1000)  GO  TO  140 
WRITE (6,640)  NGEN , KOST , NCSTF , NCNF 
WRITE<6,650) 
650  F0RMAT(15X,«CAN  NOT  PROCEEO  BECAUSE  STORAGE  IS  •  , 

1  'ALLOCATED  FOR  NO  MORE  THAN  1000  •, 

2  «CN  FEAS  LN  COM1 ) 
GC  TO  10 

C   STORE  THE  RELATIVE  COST  AND  IDENTIFYING  NUMBER 
C   FOR  THE  LATEST  CONSTRAINT  FEASIBLE  COMBINATION. 
140  CSTLC(NCNF)=COST 

NLC(NCNF)  =  2**<  Ul-iJ  +  II-1 

IF(NCNF.EQ.l)  GO  TO  160 
C   ARRANGE  VALUES  IN  CSTLC  AND  NLC  ACCORDING  TO 
C   ASCENDING  ORDER  OF  RELATIVE  COSTS. 
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NCNF1=NCNF-1 

DO    150    IL=1,NCNF1 

LI=NCNF1-IL+1 

IF(CSTLC(LH-1)  .GE.CSTLCU  I  )  )    GO    TO    160 

CHOLO=CSTLC(LI ) 

NHOLD=NLC(LI ) 

CSTLC(LI  )=CSTLC(LI*1) 

NLCtLI  )=NLC(LI-H) 

CSTLC(LI+1 )=CHOLO 
150  NLC(LI*l)=NHOLD 
160  CONTINUE 
170  CONTINUE 

IF(NCNF.EQ.O)  GO  TO  65 
C   CHECK  EACH  COMBINATION  STORED  AND  RANKED  ABOVE 
C   FOR  A  TRAVELING-SALESMAN  SOLUTION. 
180  WRITE(6,64C)  NGEN , KOST , NCSTF ,NCNF 

DC  21C  1=1, NCNF 
C   RECOMPUTE  OPTIMAL  ASSIGNMENT  BASIC-VARIABLE  VALUES. 

DO  185  K=l ,M 

XH0LD=H(K)-B(K,L1)-B(K,L2> 
185  KB(K)=XHOLD+SIGN< .5,XH0LD) 

L  =  0 

KM=M+1 

NL=NLC(I ) 
C   REGENERATE  A  COMBINATION  STORED  ABOVE. 
190  L=L+1 

LCL=M0D<NL,2) 

IF(LCL.EQ.O)  GO  TO  200 

ICL=IC(L> 
C   RECOMPUTE  NEW  BASIC-VARIABLE  VALUES 
C   FOR  LATEST  REGENERATED  COMBINATION. 

DO  195  K=1,M 

XH0LD=KB(K)-B(K,LEG(1, ICL ) ) -3( K , LEG ( 2 , ICL ) ) 
195  KB(K)=XHOLD+SIGN( .5,XH0LD) 
C   AUGMENT  ASSIGNMENT  BASIS  WITH  THE  NONBASIC 
C   VARIABLES  IN  THE  LATEST  COMBINATION. 

KM=KM+1 

IB<KM)=ICL 

KB(KM)=1 
200  NL=NL/2 

IFtNL.GT.C )  GO  TO  190 
C   CHECK  LATEST  COMBINATION  FOR  A  TRAVELING-SALESMAN  TOUR. 

COST=CSTLC(I) 

CALL  S0LCHKU220) 
C   AS  SOON  AS  A  TOUR  IS  FOUND  GO  TO  STATEMENT  220, 
C   OTHERWISE  CONTINUE. 
210  CONTINUE 

RETURN 
220  RETURN  1 

END 

SUBROUTINE  ASSIGNt*) 

COMMON  /ASGBLK/  0(1600), DB(  30>,DUAL(  80 ) , AX (  80), 
1  E, AD, W, INDEX, Kl, JS,IR,ID,ITR 

COMMON  /SHARE/  C(1600),B(  80,  80), H(  80), Z, 
1   .  LEG(2,1600) ,IB( 12^) ,MM,M,N 

WPlTE(6,35) 
35  F0RMATC1*  ) 

Kl=2 

INDEX=0 

1  =  0. 

NN  =  MM 

DO  4  1=1, M 
4  H( I)=l. 

DO  1  1=1, MM 

U=1C000. 

V=10000. 

KA=NN* I-NN+1 

KB=NN*I 

DO  6  J=KA,K8 

IF(C( J  ).LT.U)  U=C(J) 
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6      LEG(1,J)=I 

mu.EQ.O.  )  GO  TO  14 
DO  27  J=KA,KB 

27  C(J)=C(J)-U 
Z=Z-»-U*H(I) 

14     K=I+MM 

DO  3  J=ItN,NN 

IF(C( JJ.LT.V)  V=C(J) 
3      LEG(2,J)=K 

IF(V.EQ.C. )  GO  TO  1 

DO  28  J=I,N,NN 

28  C(J)=C(J)-V 
Z=Z+V*H(K) 

1  CONTINUE 
CALL  LIPU33) 
RETURN 

33  RETURN  1 
END 

SUBROUTINE  LIP{*) 

COMMON  /ASGBLK/  D<160O),DB(  30),DUAL(  8C  )  , AX <  80) 
1  EiADt Wt INDEXfKl, JStlRflDt  ITR 

COMMON  /SHARE/  C(1600),B(  80,  80), H(  80), Z, 
1  LEG(2,1600) ,IB(120) ,MM,M,N 

E=.5 

DO  5  J=1,N 

5  D(J)=-2. 
DO  3  J=1,M 
DC  2  1=1, M 

2  B(I,J)=C. 

3  B(J,J)=1. 
DO  4  I=1,M 
DB(I  )  =  -l. 
DUAL (I  )=C. 

4  IB { I )=-I 
AD=1. 
ID=1 
W=0. 

DO  6  I=1,M 

6  W=W+H(I) 
ITR=0 

7  IF(ABS(W)-E) 100,100,8 

8  DO  9  J=1,N 
IF(D(J))10,9,9 

9  CONTINUE 
GO  TO  101 

10  CX=1CC000. 
DO  11  J=1,N 
JT  =  J 

IF(DU)  )  12,11,11 

12  IF(C(J  )-.0001.LT.O. )  GO  TO  41 
IF(CX+C( J)/D(J ) ) 11 ,11, 13 

13  CX=-C( J)/D(J) 
JS=J 

11  CONTINUE 

DO  14  J=1,N 

14  C( J)=C(J)+CX*D(J) 
C( JS)=0. 
Z=Z+CX*W 

JT=JS 

DC  40  1=1, M 

40  DUALU  )=DUAL(I  )-CX*DB<  I) 

15  IF(W-E)1CC,100,41 

41  DX=0. 

DO  16  J=JT,N 

IF(C(J)-C.O0Ol  )17, 17,16 
17     C(J)=0. 

IF(D(J).GE.DX)  GO  TO  16 

DX=D(J) 
19     JS=J 

IF(D( JS) .EQ.-2.)  GO  TO  1 
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16     CONTINUE 

1      WRITE(6,213>  JS,D(JS) 

213  FORMATUH  , 3HJ S= , I  4, 3X , 2HD=, F 8.  4) 
IF(DUS) )2C,7,7 

20     CALL  ELEM 

IF( 10-1)528,26.26 
26     CALL  PIV0K&33) 

GO  TO  15 
130    CONTINUE 

INFEA  =C 

GO    TO    425 
101  WRITE     (6,2C5) 

205  FORMATUH     ,10H  INFE  AS  I  BLE  ) 

INFEA  =1 
33  RETURN  1 
425  CONTINUE 
523    1=1 

RETURN 

END 

SUBROUTINE  ELEM 

COMMON  /ASG8LK/  0(1600), DB(  80),0UAL(  80  ) , AX (  80) 
1  E,AD,W,INDEX,K1,JS,IR,ID,ITR 

COMMON  /SHARE/  C(1600),B(  80,  80), H(  80), Z, 
1  LEG(2, 1600) ,13(120) , MM, M,N 

IR=C 
J=JS 

DC  21  1=1, M 
AX( I )=0. 
DO  21  K=l,2 
K2=LEG(K,J) 
AX( I  )=AX( I  )  +  B(  I,K2) 
21  CONTINUE 
EE=E 

DO  25  1=1, M 
IF(AX(I)-E)  25,25,810 
810  IF(H( I  )-E)  82,82,23 
82  IP=I 

GO  TO  13 

23  IF(IR.EQ.O)  GO  TO  24 
IF(ABS(H(I  )-AX( I)*X)-EE)  70,70,2 

2  IF(H(I )-AX(I)*X)  24,70,25 

24  X  =  H(I)/AX(  I) 

3  IR=I 
EE=E/AX( IR ) 
GO  TO  25 

70  IF( IB( I) )  74,25,73 

73  IF(IBUR))  25,25,4 

74  IF(IB( IR) )  4,25,3 

4  IF(AX(  IR)-AX(I  ))  25,25,3 
2  5  CONTINUE 

13  RETURN 
END 

SUBROUTINE  PIVOT(*) 

COMMON  /ASGBLK/  D(1600),DB(  30 ) , DUAL (  80),AX(  80) 
1  E, AD, W, INDEX, Kl, JS,IR,ID,ITR 

COMMON  /SHARE/  C(1600),B(  80,  80), H(  80), Z, 
1  LEG(2,1600) ,IB( 120) ,MM,M,N 

84     WRITE  (6,214)  IR,AX(IR) 

214  F0RMAT(1H+,2  5X,3HIR=, I  3, 3X , 3HAX= , F 8. 4 ) 
IB(IR)=JS 

AD=AD*ABS( AX(  IR)  ) 

DI  =  IO 

ID  =  INT(DI*ABS( AX(  IR)  )*0.5) 

DI  =  ID 

IF(AX(  IR).LT.1.  +  E.AND.AX(  IR)  .GT.l.-E)  GO  TO  1 

HX=1./AX(IR) 

H( IR)=H( IR)*HX 

DO  11  J=1,M 
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11     B(  IR,J)  =  B(  IR,J)*HX 
1      DO  3  1=1, M 

IF(I.EQ.IR)  GO  TO  3 

IF(ABS(AX( I) ).LT.E)  GO  TO  3 

IF(AX(  D.LT.l.+E.AND.AXt  D.GT.l.-E)  GO  TO  4 

IF(AX(  I) .GT.-l .-E.AND.AX( D.LT.-l.+E)  GO  TO  5 

DO  2Q  J=1,M 
29     B( I.J)=B( I,J)-AX(  I )*B(  IR, J) 

H( I)=H(I)-AX(I )*H(  IR) 

GO  TO  3 

4  DO  6  J=1,M 

6  B( I,J)=B(I ,J)-B(IR,J) 
H( I)=H( I)-H( IR) 

GC  TO  3 

5  DO  7  J  =  1,M 

7  B( I. J)=B( I ,J)+B( IR,J) 
H(  I)=H(I)+H(  IR) 

3      CONTINUE 

W=W+D( JS)*H( IR) 
DO  8  J=1,M 

8  DB( J)=DB( J)-D( JS)*B< IR,J> 
E=0.5/DI 

DO  90  K=25.600,25 
IF(K-ITR)9C,91  ,95 

90  CONTINUE 
GG  TO  95 

91  DO  92  J=1,M 
DB(J)=SIGN(AINT(DI*ABS(DB(J) ) +0. 5 ) /DI , DB ( J ) ) 
DO  92  1  =  1, M 

92  B( I, J)=SIGN(AINT(DI*ABS(B( I,J))+0.5)/DI,B(I,J)) 
DO  93  1=1, M 

93  H( I)=SIGN(AINT(DI*ABS(H(I))+0.5)/DI,H( I) ) 
W=AINT(DI*W+0.5)/DI 

95     DC  34  J=1,N 
D( J)=0. 
DO  33  K=1,K1 
K2=LEG(K,J) 
IF(K2.EQ.C)  GO  TO  33 
D( J)=D(J)+DB(K2) 

33  CONTINUE 

IF(ABS(D( J))-E)35,35,34 
35     D(J)=C. 

34  CONTINUE 
ITP=ITR+1 

WRITE  (6,210)  ITRfWtZf ID, AD 
210   FORMATdH  ,5X  ,4HITR=, 15, 5X ,2HW= , F12. 8, 5X , 2HZ= , 
1    F12.5,5X.3HID=,I6,3X,3HAD=,F12.4) 
IF( ITR.EQ.300)GO  TO  21 
RETURN 
21     WRITE(6,22C) 

220    FCPMAT(//6X, 'ITERATION  LIMIT  EXCEEDED1) 
RETURN  1 
END 

SUBROUTINE  RAND(*) 

COMMON  /SHARE/  C(1600),B(  80,  80).H(  80), Z, 
1  LEG(2, 1600), IB( 120),MM,M,N 

DATA  1/3125/ 

DO  If  J=1,N 

1=^00(3125*1 ,2048) 

X=(100.*I) /2048.  +  .5 

K=X 

C( J)=K 
10  CONTINUF 

MM1=MM*1 

DO  20  J=1,N,MM1 

C( J)=999. 
20  CONTINUE 

RETURN  1 

END 
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c 
c 

c 


su 

N  = 


BROUTINE  TIMEIT(N,TIME) 

0  STARTS  CLOCK;  N=-l  STOPS  CLOCK 


10 


20 


IT 

GC 

CA 

TI 

RE 

CA 

TI 

TI 

RE 

EN 

C   PLACE 

//ASM.SY 

TIMEALL 

TIMON 


EXIT 
TIMOFF 


TOTIME 
CLOCKR 
SAVE1 
TEMPI 

C  PLACE 
//GO. FTC 
//GO. SYS 
C  PLACE 
C   PLACE 


=N+2 

TO 
LL  T 
MEM= 
TURN 
LL  T 
ME=M 
ME=( 
TURN 
D 

A  / 
SIN 

CSE 

ENT 

SAV 

USI 

LR 

ST 

LA 

L 

L 

ST 

ST 

STI 

L 

RET 

SAV 

USI 

LR 

ST 

LA 

L 

TTI 

ST 

L 

RET 

CNO 

DC 

DS 

DS 

DS 

END 

A  / 
6FCC 
IN  D 

DAT 

AN 


(2C,  10),  IT 

IMON(M) 

M 

IMOFF(M) 

TIMEM-TIME)*. 000026 


*  CARD 
DD  * 
CT 
RY 
E 


HERE 


NG 


MER 

URN 

E 

NG 


MER 


URN 
P 


TIMON, TIMOFF 

(14,12) 

TIM0N,12 

12,15 

13, TEMPI 

13,SAVE1 

2,0(1,0) 

3, TOTIME 

3, CLOCKR 

3,0(2,0) 

TASK,TUINTVL=CLOCKR 

13, TEMPI 

(14,12) ,T,RC=0 

(14,12)  ENTER 

TIMOFF, 12 

12,15 

13, TEMPI 

13, SAV El 

2,0( 1,0) 

CANCEL 

0,0(2,0) 

13, TEMPI 

(14,12),T,RC=0 

0,4 

X»7FFFFFFF» 

F 

18F 


ENTRY  VIA  -CALL  TIMON(N)-  15.5 


VIA  -CALL  TIMOFF(N)- 


HERE. 

SPACE=(CYL,( 10,1) ) 


HERE 


*  CARD 

1   DD 

D  * 

A  CARDS  HERE. 

ORANGE  /*  CARD 
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APPENDIX  B  COMPUTER  OUTPUT  FOR  A  SELECTED  PROBLEM 

2 
COST  MATRIX 

123456789  10 

999          24             18            22             31             19            33            25  30  26 

15        999            19            27            26            32            25            31  28  18 

22  23           999            23             16            29            27             18  16  27 
24          31            18          999            19            13            28              9  19  27 

23  18            34            20          999            31            24            15  25  8 


6 

24 

12 

17 

15 

10 

999 

11 

16 

21 

31 

7 

28 

15 

27 

35 

19 

18 

999 

21 

21 

19 

8 

13 

24 

18 

13 

13 

22 

25 

999 

29 

24 

9 

17 

21 

18 

24 

27 

24 

34 

31 

999 

18 

0 

18 

19 

29 

16 

23 

17 

18 

31 

23 

999 

OPTIMAL   ASSIGNMENT    SOLUTION 


BASIS 

B 

82 

0 

11 

1 

75 

1 

38 

1 

50 

1 

29 

1 

3 

0 

94 

1 

6 

1 

81 

0 

74 

0 

62 

1 

83 

1 

4 

0 

55 

0 

69 

0 

57 

1 

48 

0 

90 

0 

z 

=  140 

2 
Alternate    solution  procedures    for   this   ten-city  problem  are 

discussed    in   References   8   and    9, 
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TRY  X(  8)  =  1   REL  COST  =  0 

TRY  X(97)  =  1   REL  COST  =  1 

TRY  THE  (2**  1)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  1 
NR  WITH  REL  COST  .LE.  2  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

TRY  X(68)  =  1  REL  COST  =  2 

TRY  THE  (2**  2)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  .LE.  2  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

TRY  X(99)  =  1   REL  COST  =  2 

TRY  THE  (2**  3)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  .LE.  2   =  1 

NR  FOUND  CONSTRAINT  FEASIBLE  =     0 

TRY  X(20)  =  1   REL  COST  =  2 

TRY  THE  (2**  4)  _  i  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  . LE.  2  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

TRY  X(84)  =  1   REL  COST  =  2 

TRY  THE  (2**  5)  _  \   LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  . LE.  3  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN  TO  X(20)  AND 

TRY  THE  (2**  4)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  3  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN  TO  X(99)  AND 

TRY  THE  (2**  3)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  3  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN  TO  X(68)  AND 

TRY  THE  (2**  2)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  -  3 
NR  WITH  REL  COST  .LE.  3  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 
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RETURN  TO  X(97)  AND 

TRY  THE  (2**  1)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  1 
NR  WITH  REL  COST  ,LE.  3  =  1 
NR  FOUND  CONSTRAINT  REASIBLE  =     0 

TRY  X(58)  =  1  REL  COST  =  3 

TRY  THE  (2**  6)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  .LE.  3  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

TRY  X(  2)  =  1   REL  COST  =  3 

TRY  THE  (2**  7)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  .LE.  3  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE  =     0 

TRY  X(  9)  =  1   REL  COST  =  3 

TRY  THE  (2**  8)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  .LE.  3  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

TRY  X(19)  =  1   REL  COST  =  3 

TRY  THE  (2**  9)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  .LE.   3  =  1 

NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

TRY  X(13)  =  1  REL  COST  =  3 

TRY  THE  (2**10)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  . LE.  3  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE  =     0 

TRY  X(52)  =  1   REL  COST  =  3 

TRY  THE  (2**11)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  .LE.  3  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

TRY  X(65)  =  1   REL  COST  =  3 

TRY  THE  (2**12)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  4  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE  =     0 

RETURN  TO  X(52)  AND 

TRY  THE  (2**11)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  4  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 
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RETURN  TO  X(13)  AND 

TRY  THE  (2**10)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  4  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =      1 

RETURN  TO  X(19)  AND 

TRY  THE  (2**  9)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  4  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN  TO  X(  9)  AND 

TRY  THE  (2**  8)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  4  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN  TO  X(  2)  AND 

TRY  THE  (2**  7)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  4  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN  TO  X(58)  AND 

TRY  THE  (2**  6)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  4  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN  TO  X(84)  AND 

TRY  THE  (2**  5)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  12 
NR  WITH  REL  COST  .LE.  4  =  9 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN  TO  X(20)  AND 

TRY  THE  (2**  4)   -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  9 
NR  WITH  REL  COST  .LE.  4  =  7 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN  TO  X(99)  AND 

TRY  THE  (2**  3)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  6 
NR  WITH  REL  COST  .LE.  4  =  5 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN   TO   X(68)    AND 

TRY   THE    (2**  2)    -    1   LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  3 
NR  WITH  REL  COST  .LE.   4   =  3 

NR  FOUND  CONSTRAINT  FEASIBLE   =      0 
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RETURN  TO  X(97)  AND 

TRY  THE  (2**  1)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  1 
NR  WITH  REL  COST  .LE.  4  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

TRY  X(28)  =  1   REL  COST  =  4 

TRY  THE  (2**13)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  .LE.  4  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

TRY  X(96)  =  1  REL  COST  =  4 

TRY  THE  (2**14)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  .LE.  4  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE  =     0 

TRY  X(17)  =  1  REL  COST  =  4 

TRY  THE  (2**15)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  .LE.  4  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE  =     0 

TRY  X(92)  =  1   REL  COST  =  4 

TRY  THE  (2**16)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  5  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE  =     0 

RETURN  TO  X(17^  AND 

TRY  THE  (2**15)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  5  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN  TO  X(96)  AND 

TRY  THE  (2**14)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  5  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN  TO  X(28)  AND 

TRY  THE  (2**13)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  5  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN  TO  X(65)  AND 

TRY  THE  (2**12)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  16 
NR  WITH  REL  COST  .LE.  5  =  11 
NR  FOUND  CONSTRAINT  FEASIBLE  =     0 
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RETURN  TO  X(52)  AND 

TRY  THE  (2**11)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  16 
NR  WITH  REL  COST  .LE.   5  =  11 

NR  FOUND  CONSTRAINT  FEASIBLE   =      0 

RETURN  TO  X(13)  AND 

TRY  THE  (2**10)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  16 
NR  WITH  REL  COST  .LE.  5  =  11 
NR  FOUND  CONSTRAINT  FEASIBLE  =      1 

RETURN  TO  X(19)  AND 

TRY  THE  (2**  9)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  16 
NR  WITH  REL  COST  .LE.  5  =  11 
NR  FOUND  CONSTRAINT  FEASIBLE   =      0 

RETURN  TO  X(  9)  AND 

TRY  THE  (2**  8)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =16 
NR  WITH  REL  COST  .LE.  5  =  11 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN  TO  X(  2)  AND 

TRY  THE  (2**  7)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  16 
NR  WITH  REL  COST  .LE.  5  =  11 
NR  FOUND  CONSTRAINT  FEASIBLE   =      0 

RETURN  TO  X(58)  AND 

TRY  THE  (2**  6)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  15 
NR  WITH  REL  COST  .LE.  5  =  11 
NR  FOUND  CONSTRAINT  FEASIBLE   =      0 

RETURN  TO  X(84)  AND 

TRY  THE  (2**  5)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  17 
NR  WITH  REL  COST  .LE.  5  =  15 
NR  FOUND  CONSTRAINT  FEASIBLE   =      0 

RETURN  TO  X(20)  AND 

TRY  THE  (2**  4)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  12 
NR  WITH  REL  COST  .LE.  5  =  11 
NR  FOUND  CONSTRAINT  FEASIBLE   =      0 

RETURN  TO  X(99)  AND 

TRY  THE  (2**  3)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  7 
NR  WITH  REL  COST  .LE.   5   =  7 

NR  FOUND  CONSTRAINT  FEASIBLE   =      0 
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RETURN  TO  X(68)  AND 

TRY  THE  (2**  2)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  3 
NR  WITH  REL  COST  .LE.  5  3 

NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN  TO  X(97)  AND 

TRY  THE  (2**  1)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  1 
NR  WITH  REL  COST  .LE.  5  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

TRY  X(86)   =   1  REL  COST  =  5 

TRY  THE  (2**17)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  .LE.   5  =  1 

NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

TRY  X(54)  =  1   REL  COST  =  5 

TRY  THE  (2**18)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  .LE.  5  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

TRY  X(66)  =  1   REL  COST  =  5 

TRY  THE  (2**19)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  .LE.  5  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

TRY  X(71)  =  1   REL  COST  =  5 

TRY  THE  (2**20),  -  1  LN,  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  .LE.  5  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE  =     0 

TRY  X(25)  =  1   REL  COST  =  5 

TRY  THE  (2**21)   -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  2 
NR  WITH  REL  COST  .LE.  5  =  1 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

TRY  X(85)  =  1   REL  COST  =  5 
TRY  THE  (2**22).  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  6  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE  =     0 

RETURN  TO  X(25)  AND 

TRY  THE  (2**21)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  6  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE  =     0 
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RETURN  TO  X(71)  AND 

TRY  THE  (2**20)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  6  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =      0 

RETURN  TO  X(66)  AND 

TRY  THE  (2**19)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  6  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =     1 

RETURN  TO  X(54)  AND 

TRY  THE  (2**18)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  6  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =     1 

RETURN  TO  X(86)  AND 

TRY  THE  (2**17)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  4 
NR  WITH  REL  COST  .LE.  6  =  3 
NR  FOUND  CONSTRAINT  FEASIBLE   =     1 

RETURN  TO  X(92)  AND 

TRY  THE  (2**16)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =16 
NR  WITH  REL  COST  .LE.  6  =  11 
NR  FOUND  CONSTRAINT  FEASIBLE   =     0 

RETURN  TO  X(17)  AND 

TRY  THE  (2**15)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =16 
NR  WITH  REL  COST  .LE.  6  =  11 
NR  FOUND  CONSTRAINT  FEASIBLE  =     0 

RETURN  TO  X(96)  AND 

TRY  THE  (2**14)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  16 
NR  WITH  REL  COST  .LE.  6  =  11 
NR  FOUND  CONSTRAINT  FEASIBLE   =     1 

A  TRAVELING  SALESMAN  SOLUTION  IS  — 


X( 

2, 

1) 

=    1 

X( 

8, 

5) 

=    1 

X( 

4, 

8) 

=    1 

X( 

5, 

10) 

=    1 

x( 

3, 

9) 

=    1 

X( 

1, 

3) 

=    1 

X( 

7, 

2) 

=    1 

X( 

6, 

7) 

=    1 

x(io, 

6) 

=    1 

X( 

9, 

4) 

=    1 

Z 

=   146 
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RETURN  TO  X(28)  AND 

TRY  THE  (2**13)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  16 
NR  WITH  REL  COST  .LE.  6  =  11 
NR  FOUND  CONSTRAINT  FEASIBLE   =      0 

RETURN  TO  X(65)  AND 

TRY  THE  (2**12)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  40 
NR  WITH  REL  COST  .LE.  6  =  31 
NR  FOUND  CONSTRAINT  FEASIBLE  =      0 

RETURN  TO  X(52)  AND 

TRY  THE  (2**11).  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  37 
NR  WITH  REL  COST  .LE.  6  =  29 
NR  FOUND  CONSTRAINT  FEASIBLE   =      0 

RETURN  TO  X(13)  AND 

TRY  THE  (2**10)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  34 
NR  WITH  REL  COST  .LE.  6  =  27 
NR  FOUND  CONSTRAINT  FEASIBLE  =      1 

RETURN  TO  X(19)  AND 

TRY  THE  (2**  9)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  31 
NR  WITH  REL  COST  .LE.  6  =  25 
NR  FOUND  CONSTRAINT  FEASIBLE   =      0 

RETURN  TO  X(  9)  AND 

TRY  THE  (2**  8)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  28 
NR  WITH  REL  COST  .LE.  6  =  23 
NR  FOUND  CONSTRAINT  FEASIBLE  =      0 

RETURN  TO  X(  2)  AND 

TRY  THE  (2**  7)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  25 
NR  WITH  REL  COST  .LE.  6  =  21 
NR  FOUND  CONSTRAINT  FEASIBLE  =      0 

RETURN  TO  X(58)  AND 

TRY  THE  (2**  6)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  22 
NR  WITH  REL  COST  .LE.  6  =  19 
NR  FOUND  CONSTRAINT  FEASIBLE   =      0 

RETURN  TO  X(84)  AND 

TRY  THE  (2**  5)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  31 
NR  WITH  REL  COST  .LE.  6  =  21 
NR  FOUND  CONSTRAINT  FEASIBLE   =      0 
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RETURN  TO  X(20)  AND 

TRY  THE  (2**  4)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  15 
NR  WITH  REL  COST  .LE.  6  =  13 
NR  FOUND  CONSTRAINT  FEASIBLE   =       0 

RETURN  TO  X(99)  AND 

TRY  THE  (2**  3)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  7 
NR  WITH  REL  COST  .LE.   6  =  7 

NR  FOUND  CONSTRAINT  FEASIBLE   =       0 

RETURN  TO  X(68)  AND 

TRY  THE  (2**  2)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  3 
NR  WITH  REL  COST  .LE.   6  =  3 

NR  FOUND  CONSTRAINT  FEASIBLE   =      0 

RETURN  TO  X(97)  AND 

TRY  THE  (2**  1)  -  1  LN  COM 

NR  LN  COM  NECESSARY  TO  GENERATE  =  1 
NR  WITH  REL  COST.  LE.   6   =  1 

NR  FOUND  CONSTRAINT  FEASIBLE  =      0 

**Z  =   146  IS  OPTIMAL  ** 

ALGORITHM  TIME  =     1.5043  SECONDS 
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